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Abstract  of  Dissertation  Presented  to  the  Graduate  Council 
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Degree  of  Doctor  of  Philosophy 
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UNDER  UNIAXIAL  LOADING 

By 

Edward  Linde  Bernstein 

April,  1964 

Chairman:  Dr.  I.  K.  Ebcioglu 

Major  Department:  Engineering  Science  and  Mechanics 

The  deflections  of  a  rectangular  sandwich  panel  sub- 
ject to  a  thermal  gradient,  transverse  load,  and  combined 
edge  loading  were  calculated.   Two  opposite  edges  of  the 
panel  contained  rigid  inserts  and  were  simply  supported. 
The  other  two  edges  were  simply  supported,  clamped,  or  free, 
Small  deflection  theory  was  assumed,  as  well  as  an  ortho- 
tropic  weak  core.   The  effects  of  bending  in  the  faces  were 
also  included.   The  five  coupled,  linear  partial  differen- 
tial equations  governing  the  deflections  were  transformed 
into  two  sets  of  equations,  one  of  which  is  the  well-known 
equation  of  plane  stress.   The  remaining  set  was  simplified 
to  neglect  all  edge  loading  except  a  uniaxial  force  and  was 
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solved  by  the  use  of  Levy- type  infinite  series,  which  reduce 
the  equations  to  ordinary  differential  equations  with  con- 
stant coefficients.   The  convergence  of  these  series  was 
investigated.   A  standard  solution  of  these  equations  was 
obtained,  and  the  resulting  expressions  for  the  displace- 
ments were  written  in  closed  form. 

Values  obtained  from  the  solutions  were  compared  with 
experimental  values  reported  in  the  literature,  and  agree- 
ment was  found  to  be  good,  although  shortcomings  of  the  ex- 
perimental values  are  mentioned. 

The  variation  of  the  deflections  and  the  buckling 
loads  for  panels  with  varying  elastic  properties,  dimen- 
sions, and  edge  loading  were  calculated  on  the  University 
of  Florida  IBM  709  computer,  and  these  data  are  presented 
in  graphs.   The  values  of  the  buckling  loads  were  found  to 
be  somewhat  less  than  those  reported  in  a  less  exact  analy- 
sis.  It  was  noted  that,  in  general,  a  reversal  of  deflec- 
tions took  place  and  that  buckling  modes  were  characterized 
by  a  large  number  of  half-waves. 
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CHAPTER  I 


INTRODUCTION 


Historical  Background 

Legend  has  it  that  the  terra  "sandwich"  arose  in  1792 
when  the  Fourth  Earl  of  Sandwich  devised  a  convenient  way  to 
have  his  lunch  without  interfering  with  his  daily  card  game. 
Since  that  time  the  term  sandwich  has  come  to  mean  more  than 
merely  a  way  to  eat.   The  term  applies  to  any  object  con- 
structed like  a  sandwich,  that  is,  composed  of  three  layers, 
the  middle  one  being  of  a  different  material  than,  the  outer 
layers.   Sandwich  construction  is  a  type  of  composition  that 
has  recently  begun  to  achieve  wide  application.   The  con- 
struction is  characterized  by  two  thin  outer  layers 'known  as 
the  faces  and  a  thicker  inner  layer  known  as  the  core.   The 
faces  are  made  of  a  stiff  dense  material  and  resist  bending 
and  membrane  stresses,  while  the  core  is  made  of  a  light  ma- 
terial offering  resistance  mainly  to  transverse  shear  and 
compression. 

The  original,  and  still  most  widespread,  application 
of  sandwich  construction  is  in  air  and  space  vehicles;  its 


high  strength-to-weight  ratio  is  well  suited  for  this  pur- 
pose. One  of  the  earliest  applications  was  for  the  skin  of 
the  World  War  II  DeHaviland  Mosquito  bomber,  in  which  birch- 
wood  faces  and  a  solid  balsa  core  were  used.  At  the  pres- 
ent time  the  faces  are  frequently  made  of  steel  or  aluminum 
with  a  core  made  of  either  an  expanded  cellular  material  or 
a  corrugated  metallic  sheet. 

Sandwich  construction  is  also  found  in  curtain  wall 

2 
panels  for  office  buildings;   with  a  foam  plastic  core  it  is 

used  in  refrigerator  construction;  it  is  also  suitable  as  a 
hull  material  for  boats  when  a  balsa  core  is  used. 

Because  of  this  wide  use  of  sandwich  construction, 
especially  in  panel  form,  the  elastic  analysis  of  a  sandwich 
plate  is  a  timely  and  important  problem.   The  nature  of  the 
core  necessitates  the  consideration  of  several  effects  not 
present  in  the  usual  analyses  of  thin  plates.   The  strength 
of  corrugated  cores  will  obviously  depend  on  the  direction 
of  the  corrugations.   Cores  made  of  expanded  materials  are 
constructed  of  thin  strips  of  material  which  run  in  one  di- 
rection.  The  properties  of  either  type  of  core  are  thus 
orthotropic.   Because  of  the  relative  thickness  of  the  core, 
the  effect  of  shearing  deflection  is  important.   Since  cores 
contain  air  spaces,  and  thus  act  as  thermal  insulators, 


sandwich  construction  used  in  an  air  vehicle  will  be  subject 
to  a  temperature  difference  between  the  two  faces.   Thermal 
stresses  arise  from  this  temperature  gradient. 

The  earliest  study  of  a  sandwich  plate  was  performed 
by  Gough,  Elam,  and  DeBruyne;   they  considered  the  problem 
of  a  thin  sheet  supported  by  a  continuous  elastic  foundation 
in  order  to  study  "wrinkling, "  or  instability  of  the  faces 

characterized  by  a  large  number  of  buckling  waves.   Later, 

4  5 

Williams,  Leggett,  and  Hopkins   and  Leggett  and  Hopkins 

examined  the  Euler  type  of  buckling  of  sandwich  plates  as  a 
whole;  they  accounted  for  shearing  deformations  in  the  core 
by  using  the  "tilting  method,"  as  termed  by  March,  in  which 
a  line  originally  normal  to  the  middle  surface  tilts  during 
deformation  through  an  arbitrary  angle.  Use  of  an  approxi- 
mate strain  energy  method  determined  the  buckling  load. 
March  and  others  extended  this  method  to  include  the  effects 

of  orthotropic  cores  and  orthotropic  faces  of  differing 

7  8  9 
thicknesses  and  elastic  properties.  '  '    Libove  and  Batdorf 

developed  a  method  of  analysis  in  which  the  sandwich  plate 
is  treated  as  an  orthotropic  plate  with  the  effect  of  shear- 
ing deformation  included.     Reissner  developed  differential 
equations  for  the  finite  deflection  of  sandwich  plates. 
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Hoff  used  the  assumption  of  a  weak  core  which  resists 
only  vertical  shear  and  the  principle  of  virtual  displace- 
ments to  derive  exact  equations  for  the  bending  of  a  sand- 

12 
wich  plate.     He  was  able  to  solve  the  problem  for  a  simply 

supported  plate  by  assuming  that  rigid  inserts  were  set  into 

the  core  edges.   In  this  way  the  boundary  conditions  were 

altered  to  permit  a  solution. 

In  1951  Eringen  used  the  principle  of  virtual  work  to 

derive  extremely  general  equations  which  account  for  the  ef- 

13 
fects  of  compression  and  flexure  of  the  core.     In  1952 

Gerard,  by  an  order  of  magnitude  analysis,  examined  the 

limits  of  validity  of  the  simplifications  commonly  assumed 

14 
in  sandwich  plate  analysis.     Goodier  and  Hsu  considered 

15 
nonsinusoidal  buckling  modes.     In  1959  Chang  and  Ebcioglu 

extended  Hoff's  theory  to  include  an  orthotropic  core  and 

dissimilar  faces  and  found  the  buckling  load  for  a  simply 

supported  plate.     Recently,  Yu,  in  a  series  of  papers,  has 

included  the  effects  of  rotatory  inertia  and  shear  in  the 

faces.17'18'19'20 

Thermal  stresses  in  sandwich  plates  were  treated  by 

21 
Bijlaard,  who  used  a  homogeneous  plate  approximation. 

Chang  and  Ebcioglu  extended  their  approach  to  include  the 

22 
effect  of  a  constant  temperature  gradient,    and  Ebcioglu 


later  solved  the  problem  for  a  temperature  gradient  which 

23 
varies  arbitrarily  across  the  area  of  the  plate.     Utiliz- 
ing Hoff 's  approximation  for  the  boundary  conditions,  he 
presented  a  solution  for  the  bending  of  a  plate  simply  sup- 
ported on  all  edges.   Ebcioglu  et  al.  also  give  experimental 
results  for  the  deflection  due  to  a  thermal  gradient. 

Definition  of  the  Problem 

It  is  the  aim  of  the  present  analysis  (1)  to  extend 
Ebcioglu' s  work  in  order  to  develop  equations  which  include 
the  effects  of  bending  of  the  faces,  a  linear  temperature 
distribution  through  the  thickness,  and  various  boundary 
conditions,  (2)  to  solve  the  equations  for  the  exact  bound- 
ary conditions  rather  than  making  Hoff 's  approximations,  and 
(3)  to  present  graphs  illustrating  the  effects  of  various 
parameters  of  the  plate  on  the  deflection  and  buckling  load. 

A  Cartesian  coordinate  system  is  defined  for  a  rec- 
tangular sandwich  plate  as  shown  in  Figure  1.   The  middle 
plane  of  the  core  lies  in  the  xy-plane,  with  two  adjacent 
edges  lying  along  the  x-  and  y-axes;  the  z-axis  points  down- 
ward.  In  the  notational  system  employed  here  a  single  prime 
denotes  quantities  related  to  the  lower  face,  a  double  prime 
those  related  to  the  upper  face,  and  a  superposed  bar  those 
related  to  the  core. 
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The  plate  has  lateral  dimensions  a  and  b  along  the 
x-  and  y-directions  respectively,  and  the  thicknesses  of  the 
upper  face,  core,  and  lower  face  are  t" ,    t,  and  t'  respec- 
tively.  The  faces  may  differ  in  Young's  modulus  E,  but  have 
the  same  Poisson  ratio  /j.      The  core  is  assumed  to  be  homo- 
geneous and  orthotropic. 

The  temperature  change  T  of  each  face  is  defined  as 
the  change  from  some  constant  initial  reference  temperature 
and  is  arbitrarily  distributed  in  the  x-  and  y-directions; 
the  variation  in  the  z-direction  through  the  thickness  of 
each  face  is  given  by  the  following  equations: 

T'  =  T/S  +  T/ 

(1.1) 
T"  =  T/'%   +  T2" 

where  T/ ,  T/'#  Tz',    and  Tz"  are  functions  of  x  and  y  only; 

X,      =    z   -    s  (1.2) 


where 


s  =  s'    =  — (t  +  t' )  for  the  lower  face 


s  =  s"  =  -  — (t  +  t")  for  the  upper  face. 


Arbitrary  transverse  pressures  q'  and  q"  act  normal 
to  the  outer  surface  of  each  face,  while  on  the  edges  of 
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each  face  there  act  applied  axial  forces  0Nx*/  oNxy*,  and  0Ny* 
applied  moments  Mx*  Mx/*  and  My*  and  applied  vertical  shear 
forces  Q*   and  Qy*  where  each  of  these  quantities  is  defined 
per  unit  length  of  the  boundary.   The  subscripts  follow  the 
usual  convention  for  stress  notation.   These  forces  are  uni- 
formly distributed  in  the  z-direction  on  each  face  and  do 
not  act  on  the  core,  but  on  any  edge  where  external  loading 
is  applied  a  rigid  insert  is  set  into  the  core  for  its  pro- 
tection.  A  more  complete  description  of  the  insert  is  given 
in  the  section  on  boundary  conditions. 

The  Kirchoff  hypothesis  is  assumed  to  hold  for  the 
deflections  of  the  faces.   The  temperature  is  considered  to 
be  in  steady-state  equilibrium  and  the  effect  of  deformation 
on  temperature  is  neglected.   Account  is  taken  of  the  fact 
that  elastic  properties  vary  with  temperature  by  permitting 
each  face  to  have  a  different  value  of  E.   The  "weak  core" 
assumption  is  made,  that  is,  the  core  is  considered  to  be 
so  weak  in  its  resistance  to  bending  and  axial  stresses  that 
only  the  antiplane  stresses  Tz%  ,     Tzy ,    and  fZL    are  present. 
The  core  is  considered  incompressible  in  the  z-direction. 
It  is  assumed  that  neither  local  buckling  of  the  faces  nor 
bond  failure  between  the  faces  and  the  core  occurs. 


CHAPTER  II 
DERIVATION  OF  THE  EQUATIONS 

The  Strain  Energy  Expression 

A  macroscopic  approach  has  been  chosen  for  the  deri- 
vation of  the  differential  equations  governing  the  five  in- 
dependent components  of  displacement.   It  is  felt  that  this 
approach,  based  on  the  theorem  of  virtual  work,  more  clearly 
shows  the  assumptions  involved  in  formulating  a  mathematical 
model  which  adheres  to  physical  reality  yet  allows  a  solu- 
tion; this  approach  also  permits  a  more  straightforward 
derivation  of  the  boundary  conditions  than  would  the  alter- 
nate approach  based  on  the  microscopic  equilibrium  of  an 
element  of  the  plate. 

The  theorem  of  virtual  work  states  that  if  an  elastic 
body  is  in  equilibrium,  the  variation  of  the  strain  energy 
equals  the  virtual  work  of  external  forces  acting  on  the 
body.   The  expression  for  strain  energy  for  each  layer  of 
the  sandwich  plate  is  written  and  the  variation  taken  by 
applying  the  variational  calculus. 
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The  derivation  of  the  strain  energy  expressions  for 
the  faces  proceeds  from  the  consideration  that  they  are  ini- 
tially in  a  state  of  equilibrium  under  the  action  of  large 
membrane  stresses.   From  this  initial  state  the  faces  under- 
go additional  deformations.   The  approach  used  will  be  based 
on  that  of  Timoshenko    and  Herrmann. 

The  initial  components  of  stress  are  denoted  by  Sx , 
Sy ,    and  Sxy  .   They  are  assumed  not  to  vary  through  the  thick- 
ness of  each  face,  and  to  be  of  an  order  of  magnitude  higher 
than  the  stresses  produced  during  the  additional  deformation. 
The  following  initial  stress  resultants  are  defined: 

oN»  =  tS„ 

0Ny     =  tSy  (2.1) 

oN*y   =  tSXy 

Since  the  plate  is  initially  in  equilibrium,  the  fol- 
lowing equilibrium  equations  may  be  written  for  each  face: 


+ 


3x       d  y 

(2.2) 

^0Nxy      £0Ny 

v     +  -T— -   =    0 
0  x      dy 
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The  boundary  conditions  for  the  initial  state  are 


oN/=   0Ny 


(2.3a) 

o^xy      —     o"xy 


,NXy*  =  0N, 


on  y  =  0,  y  =  b,  and 


(2.3b) 


aNxy*  =  BN, 


on  x  =  0,  x  =  a,  where  0N,*,  0Ny*  and  „N  *  are  the  initially 
applied  loads  on  the  edges  of  each  face. 

The  additional  components  of  strain  are  assumed  small 
compared  to  unity;  it  is  further  assumed  that  the  elongation 
and  shear  terms  are  of  the  same  order  and  the  rotation  terms 
are  of  a  smaller  order.  However,  since  products  of  initial 
forces  and  rotation  terms  will  occur,  these  latter  terms 
will  be  retained.  Then  the  components  of  strain  relative  to 
the  initial  state  are 

£x  «e,  +  -^(oj/  +  u>yz)  (2.4a) 

cy  ~   ey    +    *(c^x2  +  uiz2)  (2.4b) 


2 


£z  =  ez  +  -^(<-V  +  ^*2)  (2.4c) 
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£xy  =  exy  -  —  uj,  ujy  (2.4d) 


tYz     ~    eYz    -  —  cjy  ojz  (2.4e) 


£zx  «  ezx  -  -^  o^  cj,  (2.4f) 


where  e,,  ey,  ez ,  exy ,  eyi ,  and  eiX  are  the  linear  strain 
components  and  cjXi    u>y,    and  u>z   are  the  rotation  components. 

Based  on  the  Kirchoff  hypothesis,  the  deflections  in 
the  faces  are  given  by 


A       >.  dw 

U  =  U  -  K  ^— 
ox 


v  =  v  -  S 


(2.5) 

dw 

*y 


where  u  and  v  represent  the  displacement  components  due  to 
the  stretching  of  the  middle  surface  in  each  face  from  the 
initial  equilibrium  state,  and  ?  is  given  by  equation  (1.2), 
The  rotations  u>z   are  considered  to  be  negligible  for 
plate  deformations.   Then  the  use  of  equations  (2.5)  in 
equations  (2.4)  produces  the  strain-displacement  relations: 


£*      =     I~  +  ^(T~)       "    *  A~I  (2.6a) 

ox         2   c)x  cDx 


dy  2    oy  dyi 
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Lv  =  h^   +  T*  +  ~  Is  -  2?  |V)  (2.6c) 


•*y 


2  dy    6x    dx  dy       dxdy 


£*     «  ^)*  +  ±<|^  (2.6d) 

2  ox      2  dy 


£**  =  0  (2„6e) 

Cyz   =  0  (2.6f) 

The  additional  components  of  stress  t*,  ty,  and  txy 
are  assumed  to  be  related  to  the  additional  strains  by  the 
Duhamel-Neuman  equations.   The  error  introduced  by  this 

assumption  is  of  the  order  of  the  ratio  of  the  initial 

Of* 

stress  to  the  elastic  modulus  of  the  face.     Furthermore, 

the  transverse  normal  stress  <rZz    will  be  neglected  in  ac- 
cordance with  classical  plate  theory.   Therefore,  the  fol- 

27 
lowing  equations  for  plane  stress    will  apply: 


tX   = -(<-*   +  yCcty)       -     - 


4-   _    E    if  f   \  E«T 

ty  = -(£/   +  yxtx)      - 

1  ~  /"  1  "/* 


t     =     E     f 

txy     i  +/u  L*y 


where  a.  is  the  coefficient  of  linear  thermal  expansion. 
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The  following  stress  resultants  are  defined  in  the 
faces.   The  membrane  forces  per  unit  length  are 


c*U 


N„    = 


t„d?  = 


Et 


J-^ 


1   -^[c)x       '     dy         2   <>x  2      dy 


-    N, 


Nv    = 


*/*. 


tyd? 


Et 


J-v^ 


1  -/*- 


«yy  dx        2   dy  2^  dx 


-    N 


r 
(2.7) 


N, 


Va 


t„y  d? 


Et 


J_ 


</2 


2(1    +/6) 


du   +    dy  +    dw    dw 
dy        dx        dx    dy 


where   the    thermal    force  per   unit    length   is 


KU  = 


Ecx 


T        1   -M 


Vz 


Td?   = 


J. 


*/i 


E«T2t 
1    -/" 


The  bending  moments  per  unit  length  are 


r*/a 


Mx    = 


Stxd?  =  -d 


■*/l 


b2 


w 


>x' 


/« 


h  w 
by2  J 


-    M- 


rVz 


Mv    = 


?tyd?  =  -d 


J. 


Va. 


bl 


w 


&y* 


c)  w 


-    M 


(2.8) 


r<4 


M, 


?tKyd?        =         ~D(1        -     yu) 


J. 


*/a 


S    w 
^x^y 
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where  the  thermal  moment  M.  r   is 


Mr  =  «D(1  +  /-0T, 


and  the  plate  modulus  is 


D  = 


Et" 


12(1  -  JUL*) 


expression 


The  strain  energy  VIF    in  the  faces  is  given  by  the 
27 


W,    =f 


J, 


t,ix    +    tytY    +    2t„y  £xy    -    2«T(tx    +    ty) 


"Y  "'y 


•x     ^      >-y  , 


dV 


J, 


(Sxcx     +    Sy£y    +    2Sxy£xy)dV 


(2.9) 


and  the  integration  is  extended  over  the  volume  of  the  body 
in  the  initial  state. 

The  stress  resultants  (2.1),  (2.7),  and  (2.8),  to- 
gether with  the  strains  (2.7)  and  the  temperature  distribu- 
tion (1.1),  are  substituted  into  expression  (2.9).   After 
integration  through  the  thickness,  there  results 


Wf  =  * 


(V2w)2  -  2(1  -  /A 


o    J 


b    vi    b   w    /  b   w  \2 
bx2-   by2  bxby 


-2(1  +  /u)ot.T,  v"2w>dxdy  + 
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t>     (■  a 


(  0nx  +  hi,  ) 

2 


du  +    1  ^W|; 
dx         2   e)x 


+  (CNy         +       iNy) 

2    y 


dv    +     1  j^Wj' 
dy         2    6y 


+    (oN^     +    JN,y) 


dv    ,    du       dw  _dw 
dx        by        dx   dy 


-  «t2  (Nx   +    Ny  )|dxdy 


(2.10) 


The  first  integral  in  (2.10)  represents  the  strain  energy  of 
bending,  while  the  second  represents  strain  energy  of  the 
stretching  of  the  middle  surface. 

Since  a  "weak"  core  that  is  inextensible  in  the  z- 
direction  is  assumed,  only  the  shear  stresses  Tzx  and  T zy 
exist.   If  the  rotation  terms  u>2   are  neglected,  the  ortho- 
tropic  stress-strain-displacement  relations  are 


\*     =    2G,£2X   =    G*(|£  +  ^) 

az        ox 


TZy         -       2Gy£Zy      -       Gy(^      +      ^) 


(2.11) 


where  Gx  and  Gy  are  the  ortho tropic  shear  moduli. 

It  follows  from  equations  (2.11)  and  from  the  assump- 
tion that  tz    =  0  that  the  displacements  in  the  core  are 
given  by  the  relations 
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u  =  U,  z  +  Uj 


v  =  V,  z  +  V2 


(2.12) 


where  the  functions  U,  (x,y),  U2(x,y),  V,  (x,y),  and  V2(x,y) 
are  found  from  the  conditions  of  continuity  of  displacements 
at  the  interfaces  between  the  core  and  the  faces.   In  this 
way 


U   -  u'  -  u"    t'   +    t"   dw 
'       t        2t    bx 


u,  = 


u/  +  u"  +  t/  -  t"  lw 
2         4    dx 


V,  = 


v7  -  v"  +  t'  .+ t"  & 


w 


2t    by 


v,  = 


v  +  v 


t'  -  t"  *w 

4    dy 


The  following  stress  resultants  are  defined  in  the 
core.   The  shear  forces  per  unit  length  are 


rV* 


Q*  = 


t        dx 


*/i 


•V* 


Qy  =     Tzydz  =  tGy( 


--  v'  -  v"   *  d 


w« 


By 


(2.13) 
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where 


t  2  2 


(2.14) 


The   strain   energy  Wc    of    the   core   is   of   the   form 


Wc     = 


(rzx<fZ)t    +   riydzy)dv 


(2.15) 


The  expressions  (2.12)  are  substituted  into  the 
orthotropic  stress-strain-displacement  relations  (2.11)  and 
then  into  the  strain  energy  term  (2.15).   After  integration 
through  the  core  thickness,  there  results 


Wc  =  ± 


rb  r  a 


'0    u0 


G.^+tV  +Gy(V,'+f^)2 


c>x 


dxdy 


(2.16) 


The  Equilibrium  Equations 

According  to  the  principle  of  virtual  work, 


8Vfl'  -  8  V  "  =  0 


(2.17) 


where  the  virtual  work  due  to  the  applied  loadings  is  as 
follows:  (1)  from  the  transverse  pressure 
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SV^    = 


rb 


•o     ^0 


(q"  -    q'  )Swdxdy 


(2)    from   the   edge    forces   on   each    face 


&V„    = 


(0Mx*Su+    0N     * Sv)dy  + 


(„NY*Sv  +    0N     *  5u)dx 


(3)  from  the  edge  moments  on  each  face 


rb 


*v„   - 


Mx*S(|^)dy 


M/s(-^)dx 
by 


3Mxy* 
~by~ 


r  o 


<Swdy  + 


^M,v 


x 


6wdx  -  2M„y* 


5w 


(4)  from  the  shear  forces  on  each  face 


rb 


SVa  = 


Qx*Swdy  + 


Qv*  Swdy 


In  forming  the  variations  use  is  made  of  the  formulas 


fi> 


rb 


0    Jo 


dF 

~  dxdy  = 

ox 


F(x,y) 


dy 
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2>F  ,  , 
-r—  dxdy 

ey 


F(x,y) 


dx 


and  of  integration  by  parts.   Equation  (2.17)  becomes 


|(d'  +  d")vvw  +  v*(m;  +  m;o  -  (rN/  +  tn;o-^ 

I  dxa 


0    ^0 


-  2(TN.;  +  ,N,/)£^  -  (TN/  +  rNy")-p  +  q'  -  q" 


_  ^[g,  (bu'     -    $Hl   +    tt  ^W)  +  G  (*V'     -  1^  +  tt  — )l|  5w 
L  *  dx     dx        dx2      y    dy     dy        dy2 J 


■{■ 


— r—   -  -  Gx  (- 


dx     dy 


♦  t  £)i  Su- 


— +  — +  G,  ( = +  t  —) 

dx      dy  t         dx 


S\x' 


^rN.y'  drNy'         _     /V/     -    v//  j    J)W 

dx  dy  t  dy 


6V 


dx  dy  r  t  dy  J 


dxdy 


dx  dx  ay  ox  dy 


3Mxy* 

+   tQ„ Q„*)Sw  +     (rNx    -    0n/)6u 

dy 


+     (TNxy     -    „Nx/)5v 


dy  + 


+ 


dw. 
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dM, 


(-My  +  My*)6(^)  +  (^  +  2 


dM«v     KT  dw 

— ^ +  7-N 

dx 


'"dy  +  rN<y  d 


dw 
x 


+  tQ 


dMxy* 


dx 


-  Qy*)6w  +  (rNxy  -  eN„y*)6u 


+  (rNy  -  aNy*)8v 


dx  +  2(Mxy  -  M  *) 


5w  =  0     (2.18) 


In  expression  (2.18)  the  total  membrane  forces  are  defined 
by 

TUX    =  0NX  +  Nx 


rNy  =  0Ny  +  Ny 


(2.19) 


rNxy  =  0Nxy  +  N*y 


The  integrals  extend  over  the  body  in  the  initial  state,  and 
virtual  displacements  are  taken  as  relative  to  that  state. 
Since  the  virtual  displacements  5u'  ,  Su",    Sv ' ,    8v", 
and  5w  are  arbitrary  within  the  plate,  their  coefficients 
in  equation  (2.18)  must  vanish.   In  this  manner  the  follow- 
ing equations  of  equilibrium  result.   In  the  first  of  these 
equations  the  terms  (N/  +  Nx")|~,  (N/ +  N/')^--!jf,  and 
(Nxy'  +  N*y")a  >      are  neglected.   The  terms  in  parentheses  are 
membrane  stresses  due  to  bending  of  the  plate  after  initial 
stress  and  are  negligible. 
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(D'    +    D")Vvw    -     (0N/  +     ,N/)-^7   "    2(0NK/  +    oN    ")J^L 

ox'  '  '     oxdy 


(.N/  +    oN/)|^   +    VZ(M/+    M/')     -    t 
oy  * 


7k    /du' 


r-    52, 


o,  (fa:  -  5SC  +  tt  {pf> 

ox  ox  ox' 


+    Gy(Sv 


oy  dy 


&v"  +   Vt  |fH, 


oy' 


+    q'     -    q"   =    0 


(2.20a) 


dN   '        dN,/        -      u/ 


)x 


oy 


t  ox 


=    0 


(2.20b) 


— -   +   — *-   +    Gx  ( 

ox  oy 


^     /U'    -    u"    ,     t   ow 


ox 


=    0 


(2.20c) 


5N"y/  .    aV      F    ,V    -  v"       *  dw.        n 
ox  oy  t  oy 


(2.20d) 


dx  oy  '  t  oy 


(2.20e) 


These  five  equations  in  five  unknowns  reduce  to  those  de- 

23 

rived  by  Ebcioglu    for  the  case  of  constant  temperature 

variation  through  the  plate  thickness  and  for  small 
deflections . 
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Equations  (2.20)  can  be  seen  to  contain  the  equations 
of  homogeneous  thin  plate  theory.   If  the  core  is  assumed  to 
be  rigid,  i.e.,  if  Gxt  =  Gyx    -»~»,  then  it  follows  from  equa- 
tions (2.11)  that  £xi  =  £yz    =  0  and  there  is  no  contribution 
from  the  core  to  the  total  strain  energy.   The  coefficients 
of  G„  and  Gy  in  equations  (2.20)  disappear  and,  except  for 
the  thermal  terms,  the  resulting  equations  are  the  classical 
ones  of  Saint-Venant,  describing  two  plates  constrained  to 
have  a  common  vertical  deflection  and  subjected  to  lateral 
and  edge  loads.   These  plates  correspond  in  their  dimensions 
and  properties  to  the  faces  of  the  sandwich  panel. 

Saint-Venant ° s  equations  also  result  if  the  thickness 
t  of  the  core  is  allowed  to  go  to  zero,  since  in  this  case 
the  integral  (2.15)  then  degenerates  to  zero.   In  either 
case  the  effect  of  the  core  is  only  to  constrain  the  two 
faces  to  have  a  common  vertical  deflection.   The  core  has 
no  stress-resisting  properties  of  its  own,  so  the  fact  that 
it  retains  its  thickness  for  the  first  case  is  of  no 
importance. 

It  is  noted  that  a  formal  application  of  these  limit- 
ing processes  to  equations  (2.20)  is  meaningless  and  pro- 
duces incorrect  results,  since  considerations  preceding  the 
derivation  of  this  equation  are  involved. 
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The  Boundary  Conditions 

The  boundary  conditions  arise  from  consideration  of 
the  line  integrals  occurring  in  the  variational  equation 
(2.18) o   These  line  integrals  may  be  made  to  vanish  by  set- 
ting either  the  variations  or  their  coefficients  equal  to 
zero.   The  boundary  conditions  are  linearized  in  the  same 
manner  as  equation  (2.20a).   After  use  of  equations  (2.3), 
they  are  found  to  be:  (1)  on  the  boundaries  y  =  0  and  y  =  b 
either 

u'  =  0  (2.21a) 


or 


either 


or 


either 


or 


either 


or 


Nx/  =  0  (2.21b) 


u"  =  0  (2.22a) 


N,/'  =  0  (2.22b) 


v'  =  0  (2.23a) 


N/  =  0  (2.23b) 


v"  =  0  (2.24a) 
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Ny"  =  0  (2.24b) 

either 

w  =  0  (2.25a) 


or 


^-(m/+  m;o  +  2  fiM,;  +  m,;')  +  (0n/+  ,n;)^ 

dy  ox         '  '    dy 

+  tQy  =  Qy«'  +  Qy*"+  ^-(M1C*/+  Mxy")  (2.25b) 

ox 


either 


=  0  (2.26a) 

dy 


or 


My'  +  M/=  My*'+  My*"  (2.26b) 

(2)  on  the  boundaries  x  =  0  and  x  =  a 
either 

u'  =  0  (2.27a) 


or 


either 


or 


either 


or 


N/  =  0  (2.27b) 


u"  =  0  (2.28a) 


Nx"=  0  (2.28b) 


v'  =  0  (2.29a) 
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either 


or 


either 


or 


N«/   =    0 


V"   =    0 


V  =  o 


w  =    0 


(2.29b) 


(2.30a) 


(2.30b) 


(2.31a) 


^-(M/+    M,")    +    2^(Mxy'    +    M„y")    +     (.N/+     .N,")^ 


»x 


dy 


+    tQ*    =    Q*'+    Qx*"+   7T-(Mxy"  +    M«yn 

dy 


dx 


either 


or 


>w 


=    0 


M/+    Mx"=    M„"  +    M/" 


(3)  for  the  corners  (0,0),  (a,0),  (0,b),  and  (a,b) 
either 

w  =  0 
or 


(2.31b) 


(2.32a) 


(2.32b) 


(2.33a) 


M 


»y 


=  Mxy* 


(2.33b) 


The  boundary  conditions  for  a  rectangular  sandwich 
panel  simply  supported  along  the  edges  perpendicular  to  the 
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y-axis  and  variously  supported  along  the  other  two  edges, 
containing  rigid  inserts  set  into  the  edges  of  the  core  that 
are  perpendicular  to  the  y-axis,  may  be  selected  from  the 
preceding  sets  of  boundary  conditions. 

A  simply  supported  edge  may  be  supported  either  on 
knife-edges  or  by  pins  of  negligible  cross-sectional  dimen- 
sions, about  which  the  panel  may  freely  rotate.   If  the 
latter  case  is  considered  for  a  sandwich  panel,  the  pins 
may  be  imagined  as  running  through  both  faces,  parallel  and 
arbitrarily  close  to  the  edges.   They  restrain  the  boundary 
from  displacement  perpendicular  to  itself,  but  permit  paral- 
lel motion.   The  restraint  in  the  perpendicular  direction 
will  cause  a  normal  stress  at  the  boundary. 

If,  on  the  other  hand,  the  knife-edge  type  of  simple 
support  is  envisioned,  the  edge  of  each  face  is  considered 
to  be  in  contact  with  a  wedge-shaped  support  along  the  line 
forming  the  vertex  of  the  wedge.   In  this  case  the  edge  is 
free  to  translate  in  directions  both  parallel  and  perpen- 
dicular to  itself. 

Finally,  the  effect  of  a  thin  rigid  insert  is  con- 
sidered.  The  insert  is  envisioned  as  being  set  into  the 
core  so  that  it  is  flush  with  the  edges  of  the  upper  and 
lower  faces  and  welded  to  these  faces  at  its  top  and  bottom 
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edges.   Its  effect  is  to  prevent  the  edges  from  having  dis- 
placements in  the  direction  parallel  to  themselves.   If  the 
insert  is  perfectly  insulated,  so  that  thermal  effects  cause 
no  expansion  of  it,  boundary  conditions  (2.21a)  and  (2.22a) 
or  (2.27a)  and  (2.28a)  will  apply.   If,  however,  thermal 
expansion  of  the  insert  is  allowed,  the  boundary  conditions 
must  allow  for  some  displacement  of  the  edges  parallel  to 
themselves.   In  general,  vertical  displacements  of  the  edges 
would  also  occur,  but  it  is  assumed  that  the  supports  pre- 
vent such  displacements.   The  extreme  case  occurs  when  the 
insert  offers  no  resistance  to  the  free  thermal  expansion  of 
the  edges.   Then  boundary  conditions  (2.29a)  and  (2.30a)  are 
replaced  by 


(2.34a) 


ry 

'    =    *'         T/dy 

J* 

%y 

"  =  ex" 

T/'dy 

w 

*/2 

(2.34b) 


and  conditions  (2.21a)  and  (2.22a)  are  replaced  by 


u'  =  «' 


ex 

T2'dx  (2.35a) 


■J", 
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r  * 

T/'dx  (2.35b) 


In  general,  the  displacement  along  the  edge  will  be  some- 
where between  zero  and  these  values,  depending  on  the  co- 
efficient of  thermal  expansion  of  the  material  of  the  insert. 
The  vertical  deflections  at  the  center  of  the  plate  were 
calculated  for  both  zero  deflection  and  the  extreme  values 
in  equations  (2.35);  only  a  slight  difference  resulted. 

The  problem  to  be  examined  here  is  a  plate  containing 
a  rigid  insert  along  the  two  edges  parallel  to  the  x-axis, 
which  are  simply  supported  on  knife-edges.   The  two  edges 
parallel  to  the  y-axis  may  be  simply  supported  by  knife- 
edges,  clamped,  or  free.   Furthermore,  the  only  external 
edge  loading  applied  to  the  face  edges  is  the  initial  load 
0Ny*  which  is  assumed  to  be  constant  along  the  edge.   No 
transverse  loading  acts  on  the  lower  face  (q'  =  0). 

The  boundary  conditions  to  be  used  are  as  follows. 
(1)  On  the  boundaries  y  =  0  and  y  =  b 


(2.36a) 


r  * 

u'    = 

ex' 

T2'dx 

J 

% 

■  X 

u"  = 

<x" 

T/dx 

J 

% 

(2.36b) 
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Ny'  =  0  (2.36c) 

Ny"=  0  (2.36d) 

w  =  0  (2.36e) 

My'  +  M/'=  0  (2.36f) 

where  the  right-hand  sides  of  equations  (2.35)  should  be  set 
equal  to  zero  if  the  inserts  do  not  expand  when  heated. 

(2)  On  the  edges  x  =  0  and  x  =  a,  for  the  case  of  simply 
supported  edges, 

N/  =  0  (2.37a) 

N/'=  0  (2.37b) 

N,/  =  0  (2.37c) 

NKy"  =  0  (2.37d) 

w  =  0  (2.37e) 

M/  +  Mx"=  0  (2.37f) 

For  the  case  of  clamped  edges 

u'  =  0  (2.38a) 

u"  =  0  (2.38b) 
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N„y'   =    0 


N   "=  0 


(2.38c) 
(2.38d) 


w  =  0 


dw 


=  0 


(2.38e) 
(2.38f) 


For  the  case  of  free  edges 


N„'  =  0 


(2.39a) 


Nx"=  0 


N*y'  =  0 


N  '  =  0 


(2.39b) 
(2.39c) 
(2.39d) 


(„NX'  +  *NX")~  +  ^(Mx'+  Mx")  +  2  #;(Mxy'  +  Mxy") 


+  tQx  =  0 


M/  +  My"=  0 


(2.39e) 


(2.39f) 


CHAPTER  III 
SOLUTION  OF  THE  EQUATIONS 

Transformation  into  an  Uncoupled  System 

The  solution  of  equations  (2.20)  is  difficult  because 
they  are  coupled  in  the  four  dependent  variables  u',  u",  v', 
and  v".   By  a  transformation  of  variables  two  sets  of  equa- 
tions more  amenable  to  solution  may  be  developed.   One  set 
resulting  from  this  transformation  is  composed  of  two  equa- 
tions in  two  variables.   These  equations  are  then  the  equa- 
tions of  plane  elasticity  including  the  effects  of  a  body 
force.   The  second  set  of  equations  may  be  solved  by  a 
separation  of  variables. 

To  effect  this  transformation,  new  displacement  com- 
ponents are  defined  as  follows: 


u  = 


v  = 


u  = 


V  = 


m 


1 

+ 
m 

m 

1 

+ 
1 

m 

1 

+ 
1 

m 

1  +  m 


(u'  -  u") 


(v'  -  v") 


(u"  +  mu') 


(v"  +  mv' ) 


(3.1) 
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where 


m  = 


E't' 
E"t" 


It  is  possible  to  give  a  physical  interpretation  of 
the  quantities  in  (3.1).   The  variation  of  u  and  the  varia- 
tion of  v  through  the  plate  thickness  are  given  by  straight 
lines  with  discontinuities  of  the  slope.   Figure  2  illus- 
trates this  variation  for  u.   If,  however,  it  is  imagined 
that  these  variations  are  given  by  lines  of  constant  slope, 
then  for  pure  bending  a  neutral  surface  exists,  which  is 
defined  by  the  quantities 


h' 

= 

—  A 

tt 

1    + 

m 

h" 

m 

tt 
1  +  m 

where  h'  and  h"  represent  the  distance  from  the  neutral 
surface  to  the  middle  planes  of  the  lower  and  upper  faces 
respectively.   Then  u  and  v  represent  the  displacements  of 
the  neutral  surface,  that  is,  since  the  neutral  surface  was 
defined  for  pure  bending,  u  and  v  represent  the  displace- 
ments due  purely  to  membrane  effects.   Furthermore,  u  and 
v  represent  the  displacements  of  the  middle  surface  of  the 
upper  face  relative  to  the  neutral  surface  of  the  plate. 
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t' 


*  u  - 


^ 

dw           A 

*x               / 

t'    dw 

*/ 

2     dx 

Figure  2.   Variation  of  displacement  u  through  the 
thickness  of  the  sandwich  panel 


t"/2 


T 

h" 


h' 


-*  N„" 


N» 


Jh. 


-►  N/ 


t'/2 


Figure  3.   Resolution  of  membrane  forces  into  a  force  and 
couple  system  at  the  "neutral"  surface 
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Substitution  of  the  transformations  (3.1)  into  the 
equations  for  membrane  force  resultants  (2.7)  yields  the 
corresponding  transformations  for  these  quantities: 


,NX  =  eNx'  +  0NX' 


„Ny  =  0N/  +  0Ny" 


0Nxy  =    0Nxy'  +  0Nxy' 

Nx  =  Nx'+  Nx" 

Ny  =  N/  +  Ny" 

Nxy  =  Nxy'  +  Nxy" 

Nr  =  N/  +  Nr" 


Hx  =   tfc   (N/  -  mNx") 
1  +  m 


—  A. 

Hy  =  -±£_(Ny'-  mNy") 
1  +  m 


Hxy  =  -^ — (Nxy'  -  mN  ") 
1  1  +  m  7 


HT  =  T^-tNr'  -  mNr") 
1  +  m 


where  t  has  been  defined  in  equation  (2.14). 

The  quantities  Nx  and  H„  may  be  interpreted  as  the 
statically  equivalent  resolution  of  the  membrane  forces  N„' 
and  Nx",  which  act  in  the  middle  planes  of  the  faces,  into  a 
force  and  couple  system  acting  at  the  neutral  surface  (see 
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Fig.  3) .   The  interpretation  of  the  other  quantities  follows 
in  a  similar  manner. 

Substitution  of  the  transformations  (3.1)  into  the 
equations  of  equilibrium  (2.20)  yields 


tf      d2w  k   d2w        r-,.2, 


(D      +    D    )V    w   -    0N„— - -   -    20N,y^-^ oN.t- 7   +   V    Mr 

dx'  dxoy  'dy6 


-  tk  ciia  |a 

L  m        ox 


i:   32w,         -    ,1  +   m  c>v        i:   d 


+   fct  TZ1]    +   Gr( 


w, 


>x 


+  tt  f^) 

m        ay  oyz 


-    q"  =    0 


(3.2a) 


E't' 


2m(l   -  /<*) 


d2u 


dzu 


d2i 


>x 


,+  <i-^+  u+/^ 


mt  ox  tt    ox 


(3.2b) 


E't' 


2m(l   -  /xl) 


dfv 

ay-2 


+  (i  -/<>!^+  <i+/^ 


a^u 

dxdy 


_    cy(i^S   v+    t|*)    =-^dHr 


mt 


dy  tt    dy 


(3.2c) 


dNx        dNxy 

+  — - — '- 

ox  dy 


=   0 


(3.3a) 


C>Nxy  dNy 

dx  dy 


=    0 


(3.3b) 
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Similarly,  the  boundary  conditions  become  as  follows 
(1)  On  the  boundaries  y  =  0  and  y  =  b 


Ny   =   0 


(3.4a) 


N<y  =  0 


(3.4b) 


u  = 


m 


1  +  m 


(ot'T/  -  oc"Tz")dx 


(3.5a) 


Hy    =    0 


(3.5b) 


w  =  0 


(3.5c) 


M/+  My"  =  0 


(3.5d) 


(2)  On  the  boundaries  x  =  0  and  x  =  a,  for  the  case  of  sim- 
ply supported  edges, 


N«  =  0 


(3.6a) 


N*y  =  0 


Hx  =  0 


Hxy  =  0 


w  =  0 


(3.6b) 
(3.7a) 
(3.7b) 
(3.7c) 
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My'   +       My"=        0 


(3.7d) 


For  the  case  of  clamped  edges 


u  =  0 


N„y   =  0 


u  =  0 


Hxy  =  0 


w  =  0 


>w 


=  0 


(3.8a) 
(3.8b) 
(3.9a) 
(3.9b) 
(3.9c) 
(3.9d) 


For   the   case   of    free   edges 


N*    =    0 


N„y     =    0 


H*   =    0 


tlxy    =    0 


M/    +       My"  =       0 


(3.10a) 
(3.10b) 
(3.10c) 


f  (M/  +    Mx")    +    2  -^-(M*/    +    Mx;')    +    tQ, 
ox  oy 

+  (X+  0n;')|^  =  o 

ox 


(3.10d) 
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Solution  of  the  Membrane  Force  Equations 

The  membrane  forces  in  the  initial  state  of  stress 
are  governed  by  equations  (1.2)  and  (1.3).   In  general,  a 
complete  solution  for  the  initial  state  of  stress  requires 
the  further  specification  of  the  compatibility  condition. 
However,  since  the  forces  0NX ,  0Ny ,  and  0N„y  have  been  as- 
sumed to  be  an  order  of  magnitude  higher  than  the  linear 
additional  stresses,  the  strain-displacement  relations  will 
include  nonlinear  terms.   As  a  result,  the  compatibility 
equation  becomes  nonlinear  and,  in  general,  a  solution  of 
the  state  of  stress  is  not  possible. 

If  the  initial  edge  forces  0NX*  0Ny*  and  0Nxy*  are 
constant,  then  the  initial  membrane  forces  are  constant. 
Since  0NX* =  0Nxy*  =  0,  there  results 


oNx  =  0Nxy  =  0 


(3.11) 

>Ny  =  0N/ 


and  therefore 


,NX  =  0N„y  =  0 

(3.12) 


0Ny  =  „Ny* 
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In  the  case  when  the  edges  of  the  plate  x  =  0  and 
x  =  a  are  clamped,  the  second  set  of  boundary  conditions 
(1.3)  is  replaced  by  the  conditions 


<,u  =  0 


oNxy  =  0 


The  solution  for  the  membrane  forces  requires  the  specifi- 
cation of  the  compatibility  condition.   If  it  is  assumed 
that  the  plate  is  clamped  after  the  initial  stress  is  ap- 
plied, then  equations  (3.12)  will  also  apply  to  this  case. 

The  membrane  forces  Nx ,  Ny,  and  Nxy  arising  from  the 
additional  stresses  satisfy  equations  (3.3)  and  boundary 
conditions  (3.4)  and  (3.6).   These  equations  are  identical 
to  those  just  discussed  for  the  initial  stresses.   However, 

since  these  forces  are  within  the  realm  of  linear  elasticity, 

27 
the  compatibility  equation  may  be  expressed  as 


V-V(NX  +  Ny)  +  (1  -/u)V2Nr  =  0 


28 
Now  the  classical  methods  for  plane  stress  problems  '  can 

be  used  for  the  solution  of  the  state  of  stress. 

If  NX*  =  Ny*  =  Nxy*  =  0  and  Nr  is  a  linear  function  of 

x  and  y,  where 
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N,  =  ,  1   (E'«'T2V  +  E"«"T,"t") 
'     1  -  /U 

then  the  equilibrium  equations,  compatibility  equation,  and 
boundary  conditions  are  satisfied  by 

N„  =  Ny  =  N,,y  =  0 

and  thus 

u  =  v  =  0 

Solution  of  the  Second  Set  of  Equations 

Once  the  membrane  forces  have  been  found  as  indicated 
in  the  preceding  section,  the  next  step  is  to  solve  the  set 
of  three  partial  differential  equations  (3.2),  which  are 
subject  to  the  boundary  conditions  (3.5)  and  (3.7),  (3.9), 
or  (3.10). 

First  the  system  is  transformed  into  dimensionless 

23 
units.   The  scheme  used  by  Ebcioglu    is  followed.   Dimen- 
sionless constants,  variables,  and  reference  quantities  are 
defined  as  follows: 


S  =  * 

a 


*'l 
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n2-     E't'  /^m2- 


b     1   -  ,0. 


/" 


V  =  -f(D'  +  D") 


oN/ 


J^,e 

Pe 

Pe 

r«    : 

'      t-'tGy 

g  = 

Gy 

/3  = 

b 
a 

/>  = 

1  +    m  i 
2 

oc  = 

2t 

The  temperature  gradient  is  nondimensionalized  by  the  intro- 
duction of  the  quantity 


R  ■  *  +^(ot/T/~oC//T^/)  (3.13) 

1  +  m 


The  quantity  Pe  has  been  defined  as  a  reference  quan- 
tity with  dimensions  of  force  per  unit  length;  it  represents 
the  buckling  load  of  a  simply  supported,  infinitely  long 
sandwich  strip  with  a  rigid  core,  calculated  under  the 
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assumption  that  the  faces  act  as  membranes.   P^  represents 
the  sum  of  the  buckling  loads  of  each  face  of  the  strip. 

Substitution  of  these  quantities  into  the  appropriate 
equations  and  boundary  conditions  yields  the  transformed 
equilibrium  equations: 


n 


b*u 


dzv 


2/*jP+    (1   -/*>JF  +  /3{1+M)'*l*1 


-    2gu 


,x£* 


m9t  «/*  aw        2/3mbre   ^ 
2p       31 


77" 


<H 


(3.14a) 


vr' 


d2u 


2^  +  ^(l-^^(l+^ 


2v 


t*«m  dw 


2mbr, 


2/o     di?  77" 


Br 


(3.14b) 


m 


27r2bG, 


^  + 


i    b¥ w  dvw 


_    /,5/j  la 
^9/    35 


4       as* 


*  2 


3v        amt    (r«k,«     -    1)    a^w        qmb         mb 


*1 


di," 


2G, 


2G, 


(3.14c) 


and   the   transformed  boundary   conditions,    written   in   terms   of 
displacement   components:     (1)    On   the   boundaries  >f   =    0   and 

n  =  1 


44 


i 


u  = 


mb 
/3(1  +  /<) 


R(5)d| 


% 


(3.15a) 


!*  +  ^|f  =  mbR 


(3.15b) 


w  =  0 


(3.15c) 


P*  d2 


w 


/T2  3r,z 


=  -M- 


(3.l5d) 


(2)  On  the  boundaries  5  =  0  and  £  =  1,  for  the  case  of  sim- 
ply supported  edges, 


*f?+ "£-■*>" 


(3.16a) 


(3.16b) 


w  =  0 


(3.16c) 


„2Pf       d2W 


=  -M- 


(3.l6d) 


For  the  case  of  clamped  edges 


u  =  0 


(3.17a) 


3u    .  dv    _ 
—  +  /3;— -  =  0 


(3.17b) 
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w  =    0 


*™  =    0 


di 


(3.17c) 
(3.17d) 


For   the   case   of    free   edges 


(3.18a) 


3^  <H 


(3.18b) 


(3.18c) 


/rab 


.a*3 


„  i    "  W  ,  »  x     O       W 

/3     ^T+     (2    -^>djyp 


1  +   m  s_~ 

+  tG,u 

m 


♦  f^M?  -  U M* 


(3.18d) 


It  is  seen  that  the  system  of  equations  (3.14)  is 
linear  and  nonhoraogeneous  and  is  accompanied  by  nonhomo- 

geneous  boundary  conditions.   By  the  substitution  of  auxil- 

29 
iary  functions    the  nonhomogeneous  boundary  conditions  on 

the  sides  t?   =  0  and  >j  =  1  can  be  reduced  to  homogeneous 

ones.   Then  it  is  possible  to  select  trigonometric  functions 

of  7p  that  satisfy  these  homogeneous  boundary  conditions. 

These  functions  may  be  employed  in  the  Levy  method   of 
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separation  of  variables  so  that  the  partial  differential 
equations  are  reduced  to  ordinary  ones.   It  is  necessary  to 
expand  the  nonhomogeneous  parts  of  the  boundary  conditions 
and  equations  in  a  series  of  these  trigonometric  functions; 
furthermore,  the  final  solution  is  expressed  in  such  series, 
and  the  convergence  of  these  series  is  investigated  in  order 
to  validate  their  differentiation. 

In  order  to  reduce  the  boundary  conditions  (3.15)  to 
homogeneous  form,  the  following  substitutions  are  made: 

u  «  u#«,n)  +  Jiff,  (*)  -  fi.(|)]  +  f0(4)         (3.19a) 


v  =  v<>  d,>?)  +  nz 


-  (i  -  >i)2 


g,  d) 

1         df,U)" 
"    2^       d|      J 

ge(l) 

1„     df0(|)" 

-    2^       d|      J 

(3.19b) 


w  =  w„ (S,n)  + 


>i(n  -  i) 


(n  +  Dh,^) 


+  (2  -  >?)h0(C) 


77" 


(3.19c) 


where  the  auxiliary  functions  denoted  by  f,  g  and  h  repre- 
sent the  following  nonhomogeneous  boundary  terms: 
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f*(!>  = 


/3(1  +  /O 


R(4,0)dS 


v2 


f,  (5)  = 


mb 


/3(1  +  y") 


R(f,l)d? 


J'/2 


gtf(«)  =  -±  mbR(4,0) 


g,  (I)  =  -i  xribR(5,l) 


h0(l)    =   -MTU,0) 
h,(«)  =  -MT(S,1) 

It  may  be  readily  verified  that  the  desired  reduction 
is  accomplished  by  means  of  substitutions  (3.19).   The  equa- 
tions governing  the  quantities  u0 ,  v0  ,  and  w0  are  then  es- 
tablished by  the  substitution  of  relations  (3.19)  into  equa- 
tions (3.14)  through  (3.18).   The  following  equations  result: 


rr' 


,^2u( 


d2U, 


d2v. 


2^tp"  +  {1-^^r  +  ^1+^^ 


-    2gu< 


mgt2oc/3  dw0         2/3mbre    ^r 


2/>       as 


IT' 


94 


+  J,  u,v) 


(3.20a) 
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r, 
rr 


d 2  v.  ,  d  2  v^  d  2  u  „ 

2  -T-T-  +   /32(l    -/"J-TTT-  +  /3(1   +  /xF 


dr} 


di 


H^ 


-    2vr 


— -^ =  ,       —  +    J2  (£  ,  ») 


(3.20b) 


m 


'*  2rr2bG, 


*   *' 


d*w, 


.v^L  +    2/3*   ^w*      +   ^£ 


3u 
-  />g/3 


n 


a/3   t  gm   d2w0  dv0         amt    (rek,e     -    1)     d2wc 


H*    "  p  *n 


^2 


qmb  mb       2 

2Gy         2Gy  r  J  ' 


(3.20c) 


The  boundary  conditions  on  £  =  0  and  5=1  become,  for  the 
case  of  simply  supported  edges, 


du0     bv0 


0 Tt    +  ^~dt  =  mbR  +  s'  (?'??) 


(3.21a) 


du_      dv0 


(3.21b) 


w0  =  Sj  U,t?) 


(3.21c) 


77  a   S|2 


=  -Mr  +    S^d,^) 


(3.21d) 


For  the  case  of  clamped  edges 


u0  ■  C,  (£,,!?) 


(3.22a) 
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3ua      bv„ 


(3.22b) 


w0  =  Sj  (?,^) 


(3.22c) 


3w„ 


y  =  c,(?,,) 


(3.22d) 


For  the  case  of  free  edges 


du. 


dv, 


+  S*~z —  =  mbR  +  S,  (|,^) 


l\  •     S~    d7] 


(3.23a) 


5u0     dv 


+  ^TT  =  S^U.tj) 


d^    ^  SI 


(3.23b) 


a^w0 


^i(/J  S|2 


52w 


-Mr  +  F3   (?,>?) 


(3.23c) 


P,/3 
7T2b 


^TF"  +  (2  "^dFa^ 


1  +  m  *  = 

tG,U0 


m 


b       ot, 


(i   dM^ 


b  dC 


+  Fy  (4,7?) 


(3.23d) 


In  the  above  boundary  conditions  the  functions  Jc   (£,>?), 
s(  (?,>?)/  C£-(§,7j),  and  Ft  (?,??),  where  i  =  1,  2,  3,  4,  repre- 
sent the  additional  terms  arising  from  the  substitutions 
(3.19) . 

The  boundary  conditions  on  >7  =  0  and  r\   =  1  are 


u0  =  0 


(3.24a) 
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|^+/3/a^=0  (3.24b) 

dr?        d? 


w0  =  0  (3.24c) 

d2w 

-V^r  =  °  (3.24d) 

It  is  seen  that  the  boundary  conditions  (3.24)  along 
the  sides  r\    -   0  and  rj    =    1  are  satisfied  by  terms  of  the  form 

Uo,     =   X,  (?)    sin   J7T7? 

Voj      =    Y,  (U    cos    jrr>7 

w0,     o   Zy  (?)    sin    J7T77 

Thus  the  following  solutions  are  assumed: 


u0  =  )         X,  (4)  sin  jnri  (3.25a) 

A  =  / 

v0    =     Y      YJ   (5)  cos  J77"^  (3.25b) 


oo 

w0  =  )         Zj   (?)  sin  jTT>7  (3.25c) 

a  =  / 

However,  before  attempting  an  infinite  series  approach, 
certain  mathematical  questions  must  be  examined. 
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Convergence  of  Series  Solutions 

Once  solutions  of  the  form  of  (3.25)  have  been  pro- 
posed, it  is  desirable  to  develop  some  assurance  that  the 
series  can  satisfy  the  required  equations,  in  other  words, 
that  term-by-term  differentiation  of  the  series  will  produce 
a  series  converging  to  a  proper  derivative. 

Such  considerations  must  not  be  thought  a  mere  mathe- 
matical superfluity.   Quantities  of  physical  interest  occur- 
ring in  the  deformation  of  a  plate  involve  second  and  third 
derivatives  of  the  displacements;  if  the  displacements  are 
represented  as  infinite  series,  it  may  happen  that  term-by- 
term  differentiation  does  not  provide  series  that  converge 
to  their  physical  counterparts.   For  example,  the  deflection 
y  of  a  vibrating  beam  of  length  rr,    simply  supported,  with  an 
applied  sinusoidal  bending  moment  at  one  end  is  represented 
by  the  following  equation: 


d 


Ix  +  a«  ^!z  =  o 


at'     dx* 

with  the  boundary  conditions 

y(0,t)  =  y(77,t)  =  0 

0(o.t)  -  0 


d*  v 

— — 4(77,t)  =  b  sin  ut 
<3x 
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If  a  solution  is  attempted  by  the  method  of  finite 
sine  transforms,  the  result  will  be  the  infinite  series 


y  =  $    r     alb?("1C(sin  wt   -  -^   sin   an' t) sin 

77  [_ ,    aJn      -   w  an'' 


nx 


which,  upon  term-by-term  differentiation,  will  be  found  to 
satisfy  neither  the  differential  equation  nor  the  nonhomo- 
geneous  boundary  condition;   nor  will  the  internal  bending 
moments  calculated  in  this  way  converge  to  their  true  values, 
In  order  to  be  assured  that  term-by-term  differentia- 
tion of  a  convergent  series,  particularly  a  Fourier  series, 

will  lead  to  another  convergent  series,  two  theorems  of  per- 
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tinence  may  be  noted.   The  first  is  as  follows. 


Theorem 


op 

1.   Let  there  be  given  a  series  }       f»(x) 


rimo 

whose  terms  are  dif f erentiable  in  the  interval 

J  =  [a,b]  and  which  converges  at  least  at  one  point 

go 

of  J.   If  the  series  J      fn'(x)  deduced  from  it  by 
term-by-term  differentiation  converges  uniformly  in 
J,  then  so  does  the  given  series.   Furthermore,  if 
T     f„(x)  =  F(x)  and  jT  f„'(x)  =  p(x)  ,  then 

F'(x)  =  p(x) 

It  is  therefore  necessary  first  to  perform  formally  the 
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term- by- term  differentiation  and  then  to  investigate  the 
uniform  convergence  of  the  result  in  order  to  determine  the 
validity  of  the  operation. 

The  second  theorem  relates  specifically  to  the  term- 
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by-term  differentiation  of  Fourier  series. 

Theorem  2.   If  f (x)  is  a  continuous  function  of 
period  2  rr  with  an  absolutely  integrable  derivative, 
then  the  Fourier  series  of  f  '  (x)  can  be  obtained 
from  the  Fourier  series  of  f (x)  by  term-by- term 
differentiation. 

A  corollary  of  this  theorem  applicable  to  the  inter- 

32 
val  (0,1)  is  given  as  follows. 

Theorem  2 (a) .  If  f(x)  is  continuous  and  absolutely 
integrable  on  (0,1),  then  the  Fourier  cosine  series 
JT  a„  cos  nnx  of  f (x)  on  this  interval  may  be  dif- 
ferentiated term  by  term,  while  the  Fourier  sine 
series  J"  b„  sin  nTrx  can  be  differentiated  term  by 
term  if  f (0)  =  f(l)  =  0. 

It  is  seen  that  the  Fourier  cosine  series  is  necessarily 
always  continuous  everywhere;  if  the  conditions  in  the  corol- 
lary are  met,  the  sine  series  is  also  continuous  everywhere. 
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Theorem  2  applies  to  functions  of  one  variable  only. 
However,  it  may  be  extended  to  series  of  the  form  of  equa- 
tions (3.25)  if  the  variable  \   is  treated  as  a  parameter. 

It  is  now  possible  to  proceed  with  the  justification 
of  the  assumed  solutions.   From  Theorem  2(a)  and  boundary 
condition  (3.24a)  the  series  for  u0  may  be  twice  differen- 
tiated with  respect  to  r\ .      Condition  (3.24a)  implies  that 
°  =  0  along  rj  =  0  and  r\   =  1,  and  thus  condition  (3.24b) 
may  be  replaced  by 

3  v« 

-P-   =  0  (3.26) 

The  cosine  series  (3.25b)  for  v0  can  be  differentiated  once 
with  respect  to  r\,  and  condition  (3.26)  implies  the  permis- 
sibility of  a  second  differentiation.  Finally,  conditions 
(3.24c)  and  (3.24d),  together  with  Theorem  2(a),  imply  that 
the  series  for  w0  may  be  differentiated  four  times  with  re- 
spect to  y\. 

It  is  concluded  that  it  is  permissible  to  perform 
all  partial  differentiations  with  respect  to  r\    required  to 
satisfy  the  differential  equations  and  boundary  conditions. 

The  validity  of  taking  derivatives  with  respect  to  £, 
according  to  Theorem  1,  can  be  established  only  after  the 
actual  solution  is  obtained  and  the  differentiations  in 
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question  performed.   Some  further  remarks  will  be  made  when 
this  stage  in  the  analysis  is  reached. 

Development  and  Solution  of  the  Ordinary 
Differential  Equations 

The  next  step  in  the  solution  involves  the  develop- 
ment of  ordinary  differential  equations  in  the  variable  5  . 
These  relations  arise  from  equating  coefficients  of  trigo- 
nometric terms  after  substituting  the  assumed  solutions 
(3.25)  into  equations  (3.20).   In  this  process  the  nonhomo- 
geneous  terms  of  equations  (3.20)  are  expanded  in  Fourier 
sine  or  cosine  series  on  the  interval  (0,1)  according  to 
the  following  formulas. 

(1)  The  Fourier  cosine  expansion  of  a  function  f (x)  defined 
on  0  i  x  ^  1  is 


f  (x)  =  -j-  +  )        a,-  COS  J77X 


where 


a,  =  2 


f(x)  cos  jtfxdx 


(2)  The  Fourier  sine  expansion  of  a  function  f (x)  defined 
onO^x^lis 
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where 


f(x)    =       )        by     sin   JT7x 
J  =  i 


b,     =    2        f(x)    sin   jtfxdx 
'o 


In  this  way  the  following  relations  may  be  written: 


=      \       K,j   {\)    sin   jrrri 


3R 


||  =     ^T    K2j  (I)    cos   jrr^  +  ^  K20  (5) 

;=' 

oo 

q  =       )        kjj  (£)    sin   J77"*!. 
./'  =  / 

oo 

V*Mr  =       V^    KVJ-  U)    sin  jrr^j 

Furthermore,  the  remaining  terms  on  the  right-hand 
sides  of  equations  (3.20)  are  given  by 

J/  (€,>?)  -   y~  H/y-  (€)  sin  jrrri 

j--' 
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J2  (?,7?)  =     Y    h*j  u)   cos  inn  +  2  Uzo  (5) 


Jj  (I.7})    =      /        Hj:    U)    sin   jrrri 


•""  "  I] 


3J 

In   a   similar  manner   the   nonhomogeneous    right-hand 
sides   of    the   boundary   conditions    (3.21)    through    (3.23)    on 
i,   =    0    and   §    =    1   are 

OO 

mbRU,>?)    +   S,  U,n)    =      V~     Baj   U)    sin   j  tt  >i  (3.27a) 


s2 


(%,n)    =     /        Bcy    (I)    cos    jvn  +  I   Bco   (|)  (3.27b) 


S3U,l)    =      Y     B'J  (§)    sin   JTT^  (3.27c) 

OO 

-Mr(|,^)    +    Sv(4,rj)    =      J"     B5J    (5)     sin    jTrn  (3.27d) 


y 


CfU,7)    =      ^        COJ    (5)    sin    jrrri  (3.27e) 


Cv  (§,>?)    =      Y     C3J    W    sin    J77"*?  (3.27f) 


■*  f 
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OO 

•Mr(!,T7)    +   FjtS,^)    =      V^     Fej-    (?)    sin   jtt>?  (3.27g) 


/3    dMT(g,n) 
b  di 


+    F„  (5f7i) 


JO 


(I)    sin   jrrn      (3.27h) 


As  a  result  of  applying  the  procedure  described  above, 
the  following  systems  of  ordinary  differential  equations  re- 
sult (D  denotes  the  operator  — -) : 

d£ 


f2r„/3' 


^3-  Da    -     [rc  (1   -/«)Ja    +    2g]      X,    - 


re/3(l   +/u)j 


TT 


DY; 


mgt  a/3  2/3mbr« 

DZ;     o   — =— =-    K/J      +    H; 


2/o 


n  i         »ij  "ij 


(3.28a) 


re/3{l  +  /u)  j 


77 


DX;    + 


re/32(l    -  /")       , 

^-   D1    -    (2rej2    +    2) 


77 


*2 


t  amjTT  2mbr 


Z;     =    H1-    K,;      +     H 


2/o  J  rr 


z         "2;     T    "ij 


(3.28b) 


P,  m/3v 


bGv  4 


D' 


irril 


P^mjv772        o(mt    (1   -    re'k,e)i2-TT 
2bGv  4 


•    Z, 


^   K.;    +  -*£-   K.-    +    H 


2G, 


vJj 


2G, 


'Vy  4iJ^ 


(3.28c) 
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■e        n*, 


ITT' 


/T(l    -/*)D'    -    1 


Y0    = 


mbre  2. 

nl       K20    +     ~2Hzo 


(3.28d) 


where  j  =  1,  2,  .... 

The  boundary  conditions  (3.21)  on  %    =  0  and  £  =  1 
become,  for  the  case  of  simply  supported  edges, 


/3   ^[^  -  /ujnYj  (I)    =  2Ba,  (£) 


(3.29a) 


dY, 
JTTX,  (£)  +  /?  —f-  =  B6;  (I) 

d£ 


(3.29b) 


Z,  (4)  =  Bej  (|) 


(3.29c) 


d2Z,  U) 
dS 


-  =  B9J  U) 


(3.29d) 


/3   d^    =  Bco  (4) 


(3.29e) 


where  j  =  1,  2,  .... 

The  boundary  conditions  for  the  clamped  and  free 
edges  similarly  follow. 

The  complementary  and  particular  parts  of  the  solu- 
tions to  equations  (3.28)  are  denoted  by  the  subscripts  c 
and  p  respectively  as  follows: 


Xj    —   VLcj     +  X  pj 


(3.30a) 
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Y,     =    YCJ     +    YpJ 


ZJ     ~    z<-j     +    ZPJ 


Y0    -    Yoc     +    Yop 


(3.30b) 
(3.30c) 
(3.30d) 


Since  equations  (3.28)  are  linear  with  constant  co- 
efficients, solutions  of  the  following  form  are  assumed: 


=  A,e 


=  Bje 


*;* 


\? 


=  c,e 


B0e 


*;« 


*.* 


(3.31a) 
(3.31b) 
(3.31c) 
(3.31d) 


Substitution  of  the  above  relations  into  the  homo- 
geneous form  of  equations  (3.28)  reveals  that  Xj    satisfies 
the  following  quartic  equation: 


(ree,  - 


t^t'l 


2p  ,  2P 


bGyt  a. 


+   gj2k,e 


+  j2re(g  -  1) 


2P, 


( ^ 

hGyr  eilo(   1  -yu 


-  1)  -  1 


e, 


+  j2k,e  (l  -  — £—) 

1  -/A 


(3.32) 


where 
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TT 


e,  =  (^)-  -  J 


Further, 


TT 
A„  =  ± 


/?  |/re(l  -/<) 

If  there  are  no  repeated  roots  of  equation  (3.32), 
the  complementary  solutions  are  given  by 


8 


-  L  h" e^' 


XCJ     =       ;        Aji  e*J,s  (3.33a) 

t  =  i 


t 


Ycj     =       )        B>t-  eAj'5  (3.33b) 


:  / 


8 


ZCJ     =       )        Cji    eXj^  (3.33c) 


C't 


Yoc     =   Bc/ex°*  +    B„  e"Xot  (3.33d) 


The  constants  B,,  and  Cji      are  related  to  the  constants  Aji 
by  the  relations 


Bji      =  by,  Aji 


where 
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bjt  = 


JTT_ 


\    r^ 

[2- 

-g(i 

+  /*>K 

2 

I 

^[re 

(1-A)j2 

+  2g] 

reg/3 

'(1- 

-y«)\- 

[29 

(rej2 

+  i)  -  i2 

re(l  + 

>-)]j 

(3.34) 


2/> 


cj't  ~ 


m/3gt  rtAj/ 


2re  fl'Xji 


2A-i 


77' 


-    (l-^)j2re  -  2§- 


re(l  +yu)  j/3AJt 


77 


byi 


The  constants  AJ4  are  determined  by  substituting  the 
complete  solutions  (3.30)  into  boundary  conditions  (3.29); 
eight  simultaneous  algebraic  equations  result,  and  their 
solution  may  be  written  in  the  form 


hji     = 


9/i 


(3.35) 


where  Aj  is  an  eight-by-eight  determinant 


whose  elements  are  given  as  follows  for  the  case  of  simply 
supported  edges: 


a«     =    (/^    "  /"J^A  )eA- 
aJ/L     =    (JTT  +  /3AA  bjL  )eA- 


aSL       =     CJL     e    " 


a6Z        _      CJl 


l7L 


=   A    2  c       e  •"■ 


,         2 


aS,L  +  l     -     ~CJL    e 


•A;t 


a* /•;  =   -c 


JL 


L7,<i*'     ~~  *JL    CJL     e 

s       2 
1  »,/.*;    -     ~  Ajj  C7i 


-*,-< 
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where   L  =    1,     3,     5,    7.       The   D,,      are   the   determinants    formed 
by   replacing   the   i-th   column   of  Aj    by   the   column   vector 


2Ba;  (1) 
2Boy  (0) 
Be,    (1)     - 

Bc,    (0)     - 

Bey   (1)     - 
B«;    (0)     - 

B,;   (1)     - 
B,;   (0)     - 


dXPj  (1) 
"   /?    — TZ +  /"jWY„;     (1) 

dl 

dX.,  (0) 
-Z9    TT^—  +  /"J^Ypy  (0) 


d*i 
J77XP7  (1)    -  /3 

jnXpj  (0)    -  /3 
ZPJ  (1) 

zw-  (0) 

d%;  (1) 
dp 

dzZPj  (0) 


dY^  (1) 
d4 

dYPj  (0) 
d| 


d«* 


Also, 


B 


-  A, 


Bco   (0)    -   /3 


dY, 


^1    -    B.. 

i4        J 


(i)  ♦  a  dY" (1> 


cU 


or 


/3^0(e-A°     -    eA°  ) 


Bco  (0)    -    ft 


B, 


dYc„ (0) 
dl 


-    Bco  (1)    +    /3 


dY0,  (1) 
d£ 


/3A0(e-*°     -    eXo  ) 


For   the   case  when   the   edges  £    =    0   and  %    =    1   are 
clamped,    the   elements   of   Ay   are   given  by 


aIL      =    e 


a2l    =    1 


*3L 


=    (jn-  +  /3XjLhJL  )eA- 


a„L     =   }tt  +  /3AJL  hJL 

A; 


a.SL    -    cJL  e 


a6Z      -     C/Z 


L7i 


=    A      c      e   ■"■ 


atL        ~     AjL    CJL 
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=        ~-^L 


k7,i*( 


a2,/'/     ~      1 


aV"   =    J77  +  /3A;i  b/i 


ai,Z.'f      -      "CJL     e 


'i,/*/    =     ~CJL 
V"    =    *jl  cjL   e 


where   L  =    1,    3,     5,    7.      The   D,;     are    formed   from   the   column 
vector 


C°J 

(1) 

-    X,j  (1) 

V-Qy 

(0) 

-  xw-  (0) 

Bey 

(1) 

-  jnXpj  (l) 

Bcy 

(0) 

-    jITXf.j  (0) 

Bej 

(1) 

-    Z^y  (1) 

Be, 

(0) 

-    ZPJ  (0) 

(1) 

in) 

dZ^   (1) 
dZ„,   (0) 

-   /3 


-  fi 


dYPy  (1) 

dYFj  (0) 
d4 


d£ 


For  the  case  when  the  edges  £  =  0  and  1=1  are  free, 
there  results 


*2l 


lJi 


lVi 


lSL 


l6L 


l7i 
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(/0aAAaCy4     -/O^cj,  )eA- 
/3*XjLCjL      -  Mi'Vc^ 

/3Z\LJ-    (2   -/Oj'^A 


7rJb  L 


/£ 


-7Z 


i+_m  ~- 


m 


u 


?r2b  L 


/3 


-f<t*lS.)cAA,t}e** 

A,/-     (2    -/Oj^'Ay,]    c;t     - 


i^JS    tG, 


m 


-(t'tGjc^ 


a4,i*( 
a7,i'/ 


-(/?AA    -  //JTrb^  )e"x- 
-(/?AA    -  /tjlTbjt    ) 

(in  +   /3XJL  hJL   )e"ijL 

-(/3%L2cjt    -  /*!***  cJt  )e 

Pf/3 

772b 


/52V  -    (2  -  /u)J*w*AA]  cA     - 


1  +    m   tG, 


m 


£<*  2tG„)cAAy, 


-Aj 
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P,/3 


'■1lL't 


I^[/3V"     (2   -/<>J*"ZAy4 


772b 


m 


"   ^(t2tG,)Cj,A^ 


where   L  =    1,    3,     5,    7.      The   DJt     are    formed   from  the   column 
vector 


dX-,   (1) 

2B.,-   (1)    -   /J  ^- +  /"jffY^-   (1) 

d£ 

2Boj    (0)    -  /3  dX^(°)    +  //JttY„   (0) 


d§ 


Be;    (1)    -    jrrxp;  (1)    -  ft 


B£;   (0)    -    jr7XPy    (0)    -  ft 


dY„    (1) 
d| 

dY,,   (0) 


d5 


,1   d'Zp;    (1) 


J-  Pej   (1)    -    /?*  -   "?J""    +  Mi'rr'Zpj  (1) 


d?' 


77' 


Fe;    (0)     -     ft 


,  d*Z  pj  (0)  -    t 


Fs,    (1) 


rrlh 


fl 


d| 

idJZ^(l) 
d|3 


-      (2    -/0J*TT 


z„t  dz(l) 


d5 


+  iJ_m   ^g  (1)    +  g(VtG,)dZ:-(1) 

m  b  dC 


F,;    (0)    - 


7T2h 


,3 


^ZpAO) 


PJ 
d% 


j^-   -     (2   -//)j 


.     a  dZ(0) 


/T 


d£ 


i    +   m  /s_  ,_         fi.~2.-r-    xdZ,j  (0) 

+    X  tG*XPj  (0)    +   Mt'tGj — 

m  b  d? 


Therefore,  from  equations  (3.19),  (3.25),  (3.30), 
(3.31),  and  (3.35),  the  solutions  for  the  displacement 
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components  u,  v,  and  w  are 


u(!,t)  =  nf,  (5)  +  (1  -  >,)f„  (5) 


+ 


D 


eX-{  +  XPJ   (U 


j-.l  L*l 


sin  jrrri  (3.36a) 


v(?,»|)  =  q'' 


'        2    d£ 


-  (1  -  >?)' 


2    d? 


D 


b,7-^-  eA-?  +  Y„y  (4) 


-  r      -    i  -.  t 


COS  JTDj 


B0/eAo€  +  B„  e-Ao*   +  Yop  (?) 


(3.36b) 


77 


w(?.>?)  =  — -  >7(»7  -  1) 
6P, 


(^  +  l)h,  (I)    +  (2  -  ^)h0(U 


+ 


— i c.    J'  s 


c-  ^7 


+  Z-y  (?) 


j  ■ »   i : « ; 


sin  j/r  >7      (3.  36c) 


These  equations  completely  determine  the  displacement 
components . 

Now  that  the  solutions  have  been  completely  specified, 
further  remarks  can  be  made  about  the  validity  of  forming 
partial  derivatives  of  these  quantities  with  respect  to  ^. 
Unfortunately,  because  of  the  extreme  complexity  of  the 
terms  in  each  infinite  series,  it  is  not  possible  to 
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actually  examine  the  uniform  convergence  of  the  derivatives. 
All  that  can  be  said  is  that  actual  evaluation  of  the  series 
for  w  at  the  center  of  the  plate  indicated  that  the  terms 
converged.   Uniform  convergence  was  not  shown. 

Solution  for  a  Simplified  Case 

In  the  preceding  part  of  this  analysis  the  solution 
of  a  very  general  case  of  loading  on  a  sandwich  plate  has 
been  presented.   The  remainder  of  the  analysis  will  utilize 
certain  restrictions  which  aid  in  the  tractability  of  the 
problem,  while  still  representing  a  broad  category  of 
results. 

The  following  assumptions  are  now  made: 

1.  Bending  of  the  faces  is  neglected  (Pf  =  0)  ; 

2.  The  temperature  is  constant  through  the  thickness 
of  each  face; 

3.  The  temperature  distribution  and  the  edge  loading 
are  symmetric  with  respect  to  perpendicular  bi- 
sectors of  the  edges  of  the  plate; 

4.  The  transverse  loading  q"  is  constant. 

As  a  result  of  these  restrictions,  the  following  re- 
lations hold: 

■  i 


£'(5)  -£'u)  =^TT^) 


RU,l)d!     (3.37a) 
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9/  U)  =  9o  U)  =  mbR 


h,  (?)  =  h0U)    =  0 


(3.37b) 
(3.37c) 


Mr  =  0 


(3.37d) 


The  particular  solutions  can  be  determined  from  equa- 
tions (3.28)  after  the  following  relations  are  noted: 


T,(S,7)    =    - 


re/3mb.       2  .  ,9R.t    ,, 

a      (- A  +    1   +/<)r-(S,l) 


+    2g 


mb 


/3(1   +>u) 


RU,l)cH 


Vi 


(3.38a) 


re/3 


/"       xd    R 


Titt.l)    =  "  T^rd  -/«)  (2*1  -   Dmb(l  -  -— )^-£ 

277  1    +   JU    Ot, 


+     (2>j  -    1) 


mb(l 


A< 


1   +yu 


)R(£,D 


(3.38b) 


Tj    (5,7?)      =      /> 


mb(l   +    g  1 ^— )R(£,1) 

1       +yU  1       +   /6< 


(3.38c) 


and  after  the  following  series  expansions  are  noted: 


(3.39a) 
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i    i  r  i 

2    "    TT*      /_       j« 


cos    j^r? 


y>  »,j. 


(3.39b) 


&R(S,»?) 


rl 


=     2 


3RU,»t) 


sin   jtt>7  d>7   sin   jtttj  (3.39c) 


j'i,i      "o 


07, 


cos    jn-^  dr^  cos    jTTrj  (3.39d) 


y»f,j...   Jo 


where  the  series  in  the  last  two  equations  contain  only 
terms  with  an  odd-numbered  index  because  of  the  symmetry  of 
R(£,>7).   Thus  from  equations  (3.38)  and  (3.39)  there  result 


r  i 


K,,      =    2 


N/; 


dR 

y|-  sin  jnrj  d*, 


(3.40a) 


rl 


K2J     =    2 


dR 


COS     JTTT}    6.7] 


(3.40b) 


K 


jj 


77J 


(3.40c) 


KW    =    0 


K„     =    0 


(3.40d) 
(3.40e) 
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H. 


4 

TT 


r^  iribC 


TT 


1   +  M 


dR(^l) 
-/*+    1   +  /<)      ^1 


2g 


mb 


/3(1   +  /«)    J 


I 


R(S,l)d5 


'/z 


1 
J  j 


(3.40f) 


H 


2J 


■-^(^.>»V-»-,(i.T^r)ii|£ 


/<  xd*R(£,l) 

2 


mb(l   - 


A 


1  +  // 


)R(4,1) 


(3.40g) 


hj;     =  „  fi 


mb(l  +    g 


M 


1    +  /u  1    +/" 


)RU,D 


1 
J 


(3.40h) 


H,„    =    0 


(3.40i) 


where   j    =    1,     3,    5 

Then  the  particular  solutions  XPJ  ,    Y pj  ,    and  Zpy  are 
found  from  equations  (3.28)  after  Pf    has  been  set  equal  to 
zero  and  expressions  (3.40)  have  been  substituted  for  the 
quantities  on  the  right-hand  sides  of  these  equations. 
Furthermore,  it  is  found  that  for  any  symmetrical  tempera- 
ture distribution 


Yop    =    0 
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The  particular  solutions  for  two  specific  temperature 
distributions  were  found.   For  a  constant  distribution  R  =  R„ 
these  solutions  are 

XpJ    =    Aj\    +    By 

YPy     =   C;  (3.41) 


ZPJ     =    D; 


where 


A,    = __8mbg __^ ^    ^ 


'  77/3(1   +  /<)j   [re(l   -  yu)i2    +    2g] 


B;    =    -   -i   Ay 
J  2      J 


c     =   y>  ?/   -  A  */ 


J      '   *J  ?>   "  A"  «/ 


and 
«,     =    2(rej2   +    1) 

=   t^mjjT 
2/. 
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*J  TT*)2 


mb(l — — )R 

1    +/4 


8rembgRc 


TT 


'[re(l   -//)j2   +    2g] 


&j    =    />jrr 


*J   - 


rfmtMl   -    rek,e  )7Taj* 


A/   -A 


=—   +  p 

2GV        r 


mb(l   +    ?   ~^)R0 


1   + 


Z4 


8/ambg^R 


77j(l   +/u)[re(l   -/<)ja    +    2g] 

Here   j=l,     3,    5,     .... 

For   a   double   sinusoidal   distribution   of   the    form 
R  =    R,    sin  77- ^    sin   tty\,    the   solutions    for   j    =    1    are 

x/»>    =   Ay     cos   ""  ^ 


Y,,,-    =   Ej    sin  rr{   +    F; 


(3.42) 


where 


zpJ    ~   H>    sin  77- 1    +    Ij 


A;    -    °j 


E;     =     9, 


p2; 


H, ■    =    © 


7    Py 
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Here 


where 


2mbrf 
*J  =  "  —FT- 


Vj     is    the   determinant 


*J         *j         fli 


0  P  UJ 


\)j    =    2^/J1  +    re  (1  -  //)    +    2g 


-x,  =  re/?(l  +//) 


rj     =  /og/9 


V;     =    r,/^!   -//)    +    2re   +    2 


Py     = 


mgtza/3 


2P 


TT 


77«t'm 


CJ;     = 


(,0'g  -   rek,e    +    1) 


and  P«y     is    the   determinant   formed  by   replacing   the   i-th 
column  of   Py    by   the   column  vector 

ft 

1 
0 
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Further, 


/3j 


Fj     "   "Aj«,    ^     ■  fijSj 


where 


«, 


lj         Aj  <*j    «j  -    4  5, 


a,  =  4  ^ 


7TJ       2Gy 


For   j    =    3,     5,7,     ... 


A;    =    Ej     =    H,     =    0 


The   coefficients    in    (3.27)    are   determined   from   equations 
(3.37)    as 


B0J    (0)    =    Bay    (1)    =  |  jmb[//  +     (1   -  ^)]r(1,1)      -i 


+    2mb 


R(l,77)     sin    JTT  >7  drj 


(3.43a) 


BC7    (1)    =    -Bej    (0)    =     *   ^nb  i^ll(l    -  -^-)^r 


(3.43b) 


Bej    (1)     =     Bej    (0)     =     0 


Bsy    (1)    =    B^    (0)    =    0 


B„   (1)    =    Bco  (0)    =    0 


(3.43c) 
(3.43d) 
(3.43e) 
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Caj    (1)    =    Caj    (0)    =    0 


(3.43f) 


C.J  (1)    =    C-;   (0)    =    0 


(3.43g) 


Fej    (1)    =    F„    (0)    =    0 


(3.43h) 


Fsj    (1)    =    ~F9J   (0)    = 


4(1   +    m)tGxb 
7r/3(l   +  M) 


R(§,1)(H 


1 
J 


(3.431) 


where    j=l,     3,     5,     .... 

The  assumption  of  a  temperature  distribution  sym- 
metric in  the  ^-direction  makes  possible  a  simplification 
of  the  complementary  solutions  (3.33).   The  function  Ej-(|) 
is  defined  as 


Ey-U)  -  e 


—  o  J 


+    e 


x,U-  '/*) 


and  is  seen  to  be  even  with  respect  to  the  line  £  =  1/2, 
which  bisects  the  plate;  the  function 


0;(5)  =  a*""*' 


-  e  J 


A,U-'^J 


is  seen  to  be  odd  with  respect  to  \   =   1/2.   Each  pair  of 
terms  of  the  form 


AJL  eXjii      +   A7>(„  e-Xj* 


that  occurs  in  the  complementary  solutions  (3.33)  may  be 
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written  in  terms  of  the  newly  defined  functions  as 


'it  •*, 


'/,  Aj 


'/«*. 


-  '/l  \j 


MAj,   e'"*'   +   AJjif1   e""<  )EJ    +  ±(AJt   e'"-    -   Ay>,    e  "  ^  )0 


Under  the  present  assumptions,  the  horizontal  dis- 
placement component  u  will  be  odd  with  respect  to  5  =  1/2, 
while  the  horizontal  displacement  component  v  and  the  verti- 
cal displacement  component  w  will  be  even.   The  expression 
(3.36a)  for  u  in  terms  of  the  functions  Ej    and  07  becomes 


mb 


u  = 


/3(1    +  yu) 


•5 

RU,l)d£  + 


,D>, 


Zk-1 


'/i  *,* 


D 


^*   „-'/*  *jk 


)0;, 


sin  jrrrj 


Since  u  must  be  odd  and  the  functions  EJk     and  Oy*  are  lin- 
early independent,  the  coefficients  of  E,k     must  equal  zero; 
there  result 

DJ2     =  -e*"  Dy, 


Dj,     =  -e""  D 


J  3 


D,*   =  "e  "  °Jt 
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Therefore  the  solutions  (3.36)  become 


u  = 


mb 


/9(1   +/") 


i 


R(5,l)d£   + 

'/i 


It 


e  J        -    e  ' 


A, 


+    X 


/"> 


sin   jn-^ 


(3.44a) 


v   = 


2r)  -    1 


mb(l £-_)r(£,i)1 

1   +  /*  J 


oo     r    3 

J --1,3..       i't 


byi 


D/i 


>*    +    eX-C'?'l    +    YPJ 


(3.44b) 


■     COS     j^>J 


w    = 


c,i 


Dyi 


ef*l    +    e*y«('-D 


+    Z^y    \  sin    j7T?7      (3.44c) 


/-/.J...  l<1 


Here  the  Ay,-  are  the  roots  of  the  cubic  equation 


(ree,    -   i  _  u)  (ge/  +   rek/e  j'e,    -   gj*k,e  ) 


+    j2re(g   -    1) 


eJa-  jik-(1-r^r) 


=    0 


(3.45) 


where 


e,    =   (4^)z   "  J2 


Further,  Ay  is  the  three-by- three  determinant 
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&J      =   I  ^KL 


(3.46) 


where  K,L  =  1,  2,  3;  the  elements  aKL  are,  for  the  case  when 
the  edges  \    =  0  and  E,   =  1  are  simply  supported, 

a,L     =    (eVl      +    1)  (/3AA    -  //j^   ) 

a2L    =     (eA-      -    1)  (jrr   +   /3A„  byi  ) 

a3L    =    (eA^     +    l)cA 

where  L  =  1,  2,  3.   The  D,t   are  the  determinants  formed  by 
replacing  the  i-th  column  of  Aj  by  the  column  vector 


2B0,   (1)    -  /3 


dX„;  (1) 


L/»J 


d£ 


+    /tjlTYpj     (1) 


Bcy    (1)     -    }rrXpj   (1)     -  /3 


-Z^    (1) 


Mpj   (1) 


(3.47) 


where  the  coefficients  Baj     and  BCJ     are  given  by  expressions 
(3.43a)  and  (3.43b) . 

For  the  case  when  the  edges  \    =  0  and  %   =  1  are 
clamped,  the  elements  of  A7  are 


aM  -  e""   -  1 


a2<  =  (e  "   -  1)  (jrr  +  /3A^b„  ) 
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aJt  =  (ex-  +  l)c,-t 


and  the  Dy,   are  formed  from  the  column  vector 


-Xp;  (1) 


Bcj  (1)  -  JTTXp,-  (1)  -  /3 


-ZPJ   (1) 


dY,y  (1) 


For  the  case  when  the  edges  £  =  0  and  £  =  1  are  free, 
the  elements  of  A;  are 


a„     =    (eA'L     +    1)  (/3XJL    -  /u]nhJL  ) 


a2L    =    (e  "      -    1)  (J77-  +  /3AjLbJL  ) 


a3L    =    (eA-      -    1) 


*Ji 


m  b 


and  the  DJL     are  formed  from  the  column  vector 


dX„.  (1) 
2Ba,   (1)    -  /3  =- +  /UjrrYw  (1) 


d4 


B«y  (1)    -    J7TXW  (1)    -   /3 


dY,>   (1) 


FgJ      1) tG.X^-  (1)    -  HtztG.)     — 

"J  m  b  d£ 
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The  coefficients  F9j     are  given  by  expression  (3.43i).   The 
quantities  hJl     and  c^,   are  given  by  expressions  (3.34) 

An  observation  can  be  made  on  the  roots  of  the  char- 
acteristic equation  (3.45).   If  the  equation  is  expanded 
into  a  cubic  equation  of  standard  form 

A3  +  aA1  +  bA+  c  =  0  (3.48) 

where 

A=  A/ 

it  is  seen  from  the  definitions  and  physical  meaning  of  the 
coefficients  that 

a  <  0 

b  >  0 


c  <  0    for   k,e  < 


1  +  rej2 

1  +  rej* 
From  Descartes'  rule  of  signs  it  follows  that  the  number  of 


c  >  0    for   k,P  > 


positive  real  roots  of  the  cubic  equation  will  be  either  one 

-r  or  zero  or  two  if  k,»  > — 

1  +  rej2  ,e  1  +  rej' 

In  general,  then,  as  ~kfe  ,    which  represents  the  compressive 


or  three  if  k/e  < ~r   or  zero  or  two  if  k/e  >         z  . 


edge  loading,  increases  past  this  value,  complex  values  of 
Xj    are  to  be  expected. 
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A  further  observation  is  that  when  0Nvy  =  0,  i.e., 
no  compressive  edge  loading  is  present,  a  double  root  for  Aj 
occurs;  in  this  case  the  solutions  (3.44)  approach  the  fol- 


lowing limits  as  0N 


yy 


0: 


mb 


u  = 


/9(1  +  M) 


•I 
R(S,l)dS  + 


D 


ji 


£eXj'?  -  (1  -  5)e 


\,,('-5> 


j-L> 


L't 


A,  L 


e**1  -  eXJil"" 


+   XPJ   >  sin  jirrj 


(3.49a) 


v  =  (2n  -    1)^(1  -  — ±L- )RK,1) 
2      1  +  /* 


j.i.y 


3-1*1 


5eA"?  +  (1  -  5)e 


X;,(«-J) 


-  [eX''«  +  e^""] 


2 


by/ 


**>«*  +  e*J 


J  A, 


+  Y 


pj 


cos  jf-xj 


(3.49b) 


"I 


-jf 


5ex;'   +  (1  -  §)e 


A,f('-*J 


7'/,J 
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+  9j,      e  J'       +    e 


D,-3 


ex-1  +  e*"0"" 


A,  +  Z" 


t  - ; 


sin  j7r>7       (3.49c) 


Here  the  Ay,-  are  the  roots  of  the  equation  and  are  easily 
found  to  be 


A  =  i2 


77      /2g  +    rej2  (1   -/") 
J1~~"  $   V         reg(l   -//) 


and  A;    represents    the   determinant 


A/      =      I    aKi 


whose  elements  are  given  by 


a„     =  /3 


1  +    (1  +  *Jt)eXj'      -/«jtf  e;f  +    (b„  +    ey,)eAj 


a;j     =   Jne  Jl    +  /3 


-(Ayie„  +    b„)    +       (by,   +    e„)Ay,    +   b>; 


.*;» 


ajj     =    <3ji  +     (cjt   +    <3ji) e 


A.i 


a;i     -s    (1   +    eXjL   )  {/3XJl    -  /UjnbjL  ) 


a2L     =    (e*-      -    1)  (jrr  +/9b7IA7i    ) 


i3L    =    <ex-     +    l)cJL 
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where  L  =  1,  2.   The  DJL     are  the  determinants  formed  by  re- 
placing the  i-th  column  of  A;  by  the  column  vector  (3.47). 
The  b,<   and  cJt     are  given  by  equations  (3.34); 


e;<     "  7H 


/3     rej2(3   -ju)  (1  -   g)    +    2g 


77j    reja  (1  +  /j)  (1   -   g)    -    2g 


g,<    =  - 


4/0/3 


mt^rr-'j 


^  ^  •*■ 


re  j2  (3    -  //)    +    g    rej2(l+/')    +    2 


while  XPJ  ,    Y PJ  ,    and  Z PJ     are  the  particular  solutions  of 
equations  (3.28). 


CHAPTER  IV 
COMPARISON  WITH  EXPERIMENTAL  RESULTS 

Although  there  are  many  reports  of  experiments  on 

rectangular  sandwich  panels  without  thermal  effects,  for  ex- 

33  34 

ample,  those  of  Boiler    and  Kommers  and  Norris,    the  only 

experiment  involving  thermal  deformation  known  to  the  author 

24 
is  the  one  performed  by  Ebcioglu  et  al . 

In  this  experiment  a  rectangular  sandwich  panel  was 
simply  supported  on  all  four  edges  and  a  temperature  gradi- 
ent was  applied  by  simultaneously  heating  and  cooling  oppo- 
site faces  of  the  panel.   No  compressive  load  was  applied 
and  the  deflection  was  measured  at  the  midpoint  of  the  panel, 
The  properties  of  the  panel  were  given  as  follows: 

t'  =  t"  =  0.010  in. 

t  =  0.395  in. 

/*=  0.3 

E'  =  E"  =  30xl06  psi  (4.1) 

G„  =  39,825  psi 

Gy  =  66,400  psi 

a  =  b  =  6  in. 


85 


86 


and   thus 

k/f     =    0 

Pe    =   1,853   lb/in 

t   =    1.0253 

re    =    0.067214 

oc'   =  a"  =  7.34xl0~6  in/in  °F  (4.2) 

qa  -  0 

g  =  0.6 

m  =  1 

/3  =  1 
From  Ebcioglu  et_al.  and  the  original  data  sheets  of 
the  experiment,  it  was  found  that  for  the  four  trials  made 
in  the  experiment  the  distribution  of  the  temperature  gradi- 
ent R,  defined  by  equation  (3.13),  could  be  approximated  by 
multiples  of  the  function 

-4 
R  =  (1.25  sin  nr}   +    1.373  sin  nrj  sin  rr£)10        (4.3) 

A  measure  of  the  error  caused  by  representing  the  tempera- 
ture distribution  by  equation  (4.3)  was  calculated  from  the 
formula 


E  I  a. :  -  e,-| 
Error  = (4.4) 

A  Tmax  £  i 
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where  the  aL    are  the  experimentally  recorded  values  of  the 
temperature  gradient,  the  c,;  are  the  values  of  the  gradient 
calculated  at  the  same  points  from  formula  (4.3),  and  ATmax 
is  the  maximum  difference  among  the  experimentally  recorded 
values  of  the  temperature  gradient. 

The  following  table  gives  the  temperature  gradient 
distribution,  experimentally  recorded  values  of  the  deflec- 
tion at  the  center  of  the  plate,  and  measure  of  the  error  as 
given  by  equation  (4.4)  for  each  trial. 

TABLE  1 
EXPERIMENTAL  DATA 


Trial 


Temperature 
Gradient 


Error 


Measured 
Deflection  (in.) 


R 


+0.005 
-0.056 


0.013 


1.1  R 


+0.023 
-0.089 


0.0135 


1.2  R 


+0.043 
-0.079 


0.0195 


1.8  R 


+0.065 
-0.080 


0.021 
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Since  no  edge  loading  was  applied,  the  analytical 
solution  to  this  problem  is  given  by  equation  (3.49c).   The 
particular  solutions  for  the  temperature  distribution  R  as 
given  in  equation  (4.3)  are 


Xp    =  -(0.262  cosT7?)10"3 


Yp    =  -(0.477  +  0.262  sin7r£)10-3  (4.5) 


ZP   m    (9.011  +  5.066  sin  77£)10  3 


for  j  =  1,  and 


XPJ     =  0 


Y„-  =  0 


ZPJ     =  0 


for  j  =  3,  5,  .... 

Substitution  of  the  values  given  in  (4.1)  and  (4.2) 
and  the  particular  solutions  specified  above  into  equation 
(3.49c)  produces  the  theoretical  results  given  in  Table  2 
for  the  four  trials  reported  in  the  experiment.   The  com- 
parison between  the  theoretical  and  experimental  results  is 
illustrated  in  Figure  4. 
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0.030  i- 


0.025  - 


0.000 


1.0 


1.2        1.4        1.6 

Temperature  gradient 


1.8 


2.0 


Figure  4.   Comparison  of  theoretical  and  experimental 

values  of  the  deflection  at  the 

center  of  the  plate 
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TABLE  2 
THEORETICAL  RESULTS 


Theoretically 
Trial  Calculated 

Deflection  (in.) 


% 
Error 


1  0.01335  0.0 

2  0.01468  8.9 

3  0.01602  17.9 

4  0.02403  14.3 

The  agreement  with  the  experiment  is  seen  to  be  ex- 
cellent for  trial  1,  good  for  trial  2,  and  poor  for  the  last 
two  trials.   Nevertheless,  the  results,  while  somewhat  in- 
conclusive, do  tend  to  indicate  satisfactory  agreement  with 
the  theory,  since,  as  Table  1  shows,  the  assumed  temperature 
distribution  most  closely  fits  the  actual  distribution  in 
trial  1.   Furthermore,  the  experimental  data,  as  presented 
in  the  appendix,  show  that  for  all  four  cases  there  is  some 
doubt  whether  the  experimental  temperature  distribution  was 
actually  symmetrical  in  the  y-direction;  since  data  were 
taken  at  only  two  points  in  the  half  of  the  plate  necessary 
to  establish  this  fact,  no  definite  conclusion  can  be 
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reached.   For  trials  2,  3,  and  4  the  data  show  a  definite 
lack  of  symmetry  in  the  5-direction.   It  is  felt,  therefore, 
that  the  data  of  trial  1  are  the  most  significant,  and  the 
large  errors  for  the  other  three  trials  are  not  alarming. 
The  need  for  further  experimental  evidence  is  indicated. 


CHAPTER  V 

COMPUTATION  AND  DISCUSSION  OF  THE  SOLUTION 

Computation  of  the  Solution  for  Various 
Parameters  of  the  Panel 

In  order  to  investigate  the  dependence  of  the  deflec- 
tion and  buckling  of  the  panel  on  its  properties  and  dimen- 
sions, equation  (3.44c)  was  evaluated  on  the  University  of 
Florida  IBM  709  computer  for  a  variety  of  cases.   It  was 
decided  to  limit  consideration  to  panels  with  faces  of  equal 
thickness,  which  are  simply  supported  on  all  four  edges  and 
are  not  subject  to  a  transverse  load;  in  this  way  the  number 
of  variables  was  reduced  to  manageable  size,  while  the 
ability  was  retained  to  vary  such  quantities  as  the  ratio  of 
face  thickness  to  core  thickness,  aspect  ratio  of  the  plate, 
shear  and  Young's  moduli,  and  the  temperature  distribution 
and  edge  compression. 

The  properties  and  dimensions  of  sandwich  panels  used 
in  practice  are  as  varied  as  their  applications.   Reference 
to  industrial  and  design  journals  produced  values  of  the  di- 
mensions ranging  from  t'  =  t"  =  0.006  in.  and  t  =  0.44  in., 
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as  used  in  the  wing  of  a  supersonic  bomber,    to  t  =  t"  = 
0.064  in.  and  t  =  18  in.,  as  used  in  curtain  wall  construc- 
tion.  Honeycomb  core  shear  moduli  varied  from  G*  =  9,600  psi 
and  Gy  =  21,600  psi  to  G,  =  145,000  psi  and  Gf    =  81,000  psi.36 
The  usual  facing  materials  were  steel,  with  a  Young's  modu- 
lus E  =  30x10  psi,  and  aluminum,  with  E  =  10x10  psi. 

A  useful  range  of  variation  of  the  dimensionless  pa- 
rameter re,  given  by  the  relation 

rr'E't't 


2b*(l  -  /0G, 


can  be  chosen  from  the  above  information  as 

0.01  <  re  <  2.00 

The  values  re  =  0.017  and  g  =  0.6  given  in  the  experiments 
reported  by  Ebcioglu   were  chosen  as  typical  values.   A 
typical  value  of  the  aspect  ratio  /3   was  selected  as  2.0. 
The  nondimensionalization  of  the  deflection  w  was 
accomplished  by  the  definition  of  w  as 


■»      w 
w  = 


R0b/2ta 

where  R0  is  the  maximum  value  over  the  surface  of  the  plate 
of  the  temperature  gradient  function  R  given  in  equation 
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(3.13).   The  quantity  w  was  calculated  as  a  function  of  the 
applied  edge  loading  parameter  k,e  .   The  calculations  were 
performed  for  the  two  temperature  distributions 

R  =  R0  =  constant 

R  =  R0  sin  n  4  sin  Trr\ 

for  varying  re  and  for  varying  ft . 

The  effect  on  the  deflection  of  the  sandwich  panel 
due  to  the  degree  of  lengthwise  thermal  expansion  of  the  in- 
sert was  also  examined.   Boundary  condition  (3.15a)  was  de- 
rived under  the  assumption  that  the  thermal  expansion  of  the 
edges  r\   =   0  and  rj   =  1  was  not  inhibited  by  the  insert.   At 
the  opposite  extreme  is  the  assumption  that  no  thermal  ex- 
pansion at  all  is  allowed.   The  deflection  was  calculated 
for  both  extremes  for  varying  re. 

A  computer  program,  written  in  the  FORTRAN  symbolic 
language,  was  designed  in  order  to  accomplish  the  tasks  de- 
lineated above.   This  program  instructed  the  computer  to 
evaluate  equation  (3.44c),  with  provisions  for  varying  re , 
/3  ,  g,  yu,    and  ~kre  .      An  option  was  included  so  that  those 
terms  in  the  quantities  By,  which  arise  from  the  right-hand 
side  of  boundary  condition  (3.15a)  could  be  set  equal  to 
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zero.   Furthermore,  either  of  the  particular  solutions  (3.41) 
or  (3.42)  could  be  substituted  into  the  equation,  depending 
on  the  temperature  distribution. 

The  cubic  equation  (3.45)  was  solved  by  iteration  of 
Newton's  method.   Three  forms  of  equation  (3.44c)  were  used 
in  the  program:  one  as  given;  one  with  real  and  imaginary 
parts  written  for  the  case  when  two  of  the  roots  Aj  of  equa- 
tion (3.48)  are  complex  conjugates;  and  one  with  real  and 
imaginary  parts  for  the  case  when  one  of  the  roots  Aj  is 
negative  and  Xj   =  yAj  is  thus  imaginary. 

The  infinite  series  in  the  solution  was  truncated  at 
N  terms  either  when  the  series  was  monotonically  decreasing 
and  the  N-th  term  of  the  series  was  a  given  percentage  of 
the  first  term  or  when  N  equaled  a  given  maximum  number. 

A  few  limitations  were  imposed  on  the  calculations. 
When  /3  was  less  than  0.75,  the  roots  of  the  cubic  equation 
were  so  large  that  terms  of  the  type  e  ■"  involved  in  calcu- 
lation of  the  quantity  Ay  exceeded  the  maximum  magnitude 
allowable  in  FORTRAN.   For  the  same  reason,  no  more  than 
fourteen  terms  of  the  infinite  series  could  be  calculated. 

About  fifteen  minutes  of  actual  running  time  on  the 
IBM  709  was  required  for  all  of  the  results  obtained. 
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Discussion  of  the  Results 

The  results  of  the  calculations  discussed  in  the  pre- 
ceding section  are  shown  in  Figures  5  through  8.   Each  of 
these  graphs  shows  the  dimensionless  deflection  parameter  w, 
evaluated  at  the  center  of  the  plate,  as  a  function  of  the 
edge  loading  parameter  k/e  .   In  Figure  5  the  surface  of  the 
plate  is  subject  to  a  double  sinusoidal  temperature  gradient 
distribution  given  by 

R  =  R0  sin  rrt,    sin  tt*] 

and  the  effect  of  varying  the  parameter  re    is  shown.   In 
Figure  6  the  same  information  is  presented  for  a  constant 
temperature  distribution.   Figures  7  and  8  illustrate  the 
effect  of  a  varying  aspect  ratio /J  on  a  plate  subject  to  the 
two  temperature  gradients. 

A  basic  difference  between  the  solution  for  a  double 
sinusoidal  temperature  gradient  distribution  and  for  a  con- 
stant distribution  (when  there  is  no  transverse  load  q)  is 
that  for  the  former  case  the  infinite  series  of  the  solution 
(3.44c)  reduces  to  its  first  term.   This  fact  is  seen  from 
the  particular  solutions  (3.42),  which  are  zero  for  j  >  1. 

From  Figures  5  and  6  it  is  seen  that  the  deflection 
increases  throughout  most  of  the  loading  process  as  the 
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parameter  re  increases,  while  the  buckling  load  decreases. 
A  marked  difference  between  the  two  graphs  is  observed.   The 
deflection  monotonically  increases  until  it  becomes  unbounded 
and  positive  at  the  buckling  load  in  the  case  of  a  double 
sinusoidal  temperature  gradient  distribution;  when  the  dis- 
tribution is  constant,  the  deflection  curve  starts  to  in- 
crease for  lower  values  of  k,e  ,  but  then  turns  downward  and 
decreases  until  the  deflection  becomes  negatively  unbounded 
at  a  lower  buckling  load  than  for  the  corresponding  case  in 
Figure  5.   This  effect  is  known  as  "reversal  of  deflections." 

The  difference  in  the  response  of  the  panel  to  the 
two  temperature  distributions  is  due  to  the  presence  of  all 
the  terms  in  the  infinite  series  for  the  second  case,  as 
mentioned  above.   Buckling  occurs  at  the  lowest  value  of  k/e 
so  that  the  denominator  A_,  of  a  term  in  the  series  approaches 
zero.   The  lowest  such  value  of  k,e  need  not  cause  the  de- 
nominator in  the  first  term  of  the  series  to  vanish;  the 
buckling  load  may  be  determined  by  any  term  in  the  series, 

depending  on  the  properties  of  the  plate.   This  phenomenon 

25 
occurs  in  homogeneous  thin  plate  theory  also.     In  order  to 

describe  this  effect,  the  term  "n-th  mode  buckling"  is  in- 
troduced, meaning  that  the  critical  load  occurs  when  the 
denominator  of  the  series  term  with  index  j  =  n  vanishes. 
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From  another  point  of  view,  if  the  external  forces  acting  on 
the  plate  were  due  only  to  the  compressive  edge  loading  and 
no  temperature  gradient  were  applied,  then  the  problem  would 
become  the  associated  eigenvalue  problem  of  the  one  studied 
here;  in  this  case  the  critical  load  would  cause  the  plate 
to  buckle  into  the  n-th  mode  shape,  or  into  n  half-waves  in 
the  £ -direction.   When  n  >  1,  the  effect  is  denoted  as 
"higher  mode  buckling."   The  buckling  load  is  always  of  the 
higher  mode  type  when  /3  =    2.0,  g  =  0.6,  and  re  varies  from 
0.01  through  2.00,  as  seen  in  Figures  5  and  6. 

In  the  evaluation  of  equation  (3.44c)  on  the  computer 
it  was  found  invariably  that  each  term  in  the  series  alter- 
nated in  sign.   This  result  is  not  unexpected  because  of  the 
presence  in  the  solution  of  the  factor  sin  JTrrj,    where  r[    -   1/2 
at  the  center  of  the  plate;  however,  the  sign  of  the  remain- 
ing part  of  the  solution  could  not  be  determined  beforehand 
for  the  general  case.   It  is  interesting  to  observe  in  Fig- 
ure 6  that  whenever  higher  mode  buckling  occurred,  the  mode 
involved  had  an  index  of  j  =  3,  7,  11,  ...,  so  that  the  term 
had  a  negative  sign  and  thus  the  deflection  always  became 
negatively  unbounded. 

The  deflection  at  the  center  of  the  plate  varies 
similarly  when  re  is  held  constant  and  /3  is  varied.   As  /3 
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increases,  the  deflection  for  a  given  value  of  k/e  decreases 
and  the  buckling  load  increases,  as  is  shown  in  Figures  7 
and  8.   From  Figure  8  it  is  seen  that  the  values  of  the  pa- 
rameters chosen  permit  first-mode  buckling  for  a  range  of 
values  of  ,/3.     Once  again  the  heuristic  result  is  noted  that 
higher  mode  buckling  was  characterized  by  a  negative  value 
of  deflection. 

One  of  the  results  of  the  present  analysis  that 
should  be  useful  is  the  calculation  of  more  accurate  values 
of  the  buckling  load.   When  specific  values  of  the  proper- 
ties and  dimensions  are  known,  the  buckling  load  can  be 
found  by  evaluation  on  the  computer.   Unfortunately,  the 
expression  (3.46)  for  Aj   is  so  cumbersome  that  an  analytical 
solution  for  the  buckling  load  is  forbiddingly  difficult. 

However,  the  deflections  and  buckling  loads  were  found  to 

23 


agree   closely  with   Ebcioglu's    results 
the  buckling    load   is 


His   expression   for 


k,e    =      re  (1  -  /u) 


/3b9  +  /3*r  (1  +    2g)    +  /3\T(2  +    g)    +    j 


'] 


+    2g(/?2  +    j')1     /j 


re'(l  -  /j)  (j2   +    /3l) 


(5.1) 


+    re  \/3l 


2  +   g(l  -  /u) 


+   j' 


2g  +    (1  -  //) 


+    2g 
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A  comparison  of  the  deflection  calculated  from  Eb- 
cioglu c s  results  with  the  deflection  calculated  from  the 
present  results  is  shown  in  Figure  9.   Table  3  shows  the 
comparison  of  the  buckling  loads.   The  error  is  never  great- 
er than  10  per  cent.   Therefore,  it  is  felt  that  using  equa- 
tion (5.1)  for  approximate  calculations  of  the  buckling  load 
and  its  behavior  as  a  function  of  re  ,  /3 ,    and  j  is  justified. 

TABLE  3 
COMPARISON  OF  BUCKLING  LOAD  VALUES 


re  /3 


~kfe     at 

Buckling 

Ebcioglu' s 

Present 

Results 

Results 

13.5 

12.9 

10.6 

10.1 

7.41 

7.1 

3.85 

3.9 

2.00 

2.0 

1.00 

1.0 

0.50 

0.5 

2.36 

2.4 

3.83 

3.6 

14.8 

14.1 

25.6 

23.8 

43.5 

42.7 

58.9 

59.4 

0.025  2.0 

0.050  2.0 

0.100'  2.0 

0.250  2.0 

0.500  2.0 

1.000  2.0 

2.000  2.0 

0.017  0.75 

0.017  1.00 

0.017  2.00 

0.017  3.00 

0.017  5.00 

0.017  10.00 
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Figure  9.   Comparison  of  the  deflection  at  the  center  of  a 
plate  under  a  temperature  gradient 
R  =  R0  sin  rrl   sin  Trr\   as  a  function 
of  edge  loading 
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The  buckling  mode  is  determined  by  the  equation 
dk,e 


dj 


=  0  (5.2) 


where  k/e  is  given  by  expression  (5.1).   This  equation  lo- 
cates the  value  of  j  at  which  the  smallest  critical  value  of 
k,e  occurs.   Since  j  must  be  an  integer,  this  equation  may 
only  produce  an  approximation  of  the  exact  value  of  j . 

It  is  interesting  to  note  that  cases  exist  where 
there  are  no  real  roots  of  equation  (5.2).   When  (5.2)  is 
written  in  full,  a  quartic  equation  results.   If,  for  exam- 
ple, /3  =    2.0,  g  =  0.6,  and  re  >  0.495,  all  coefficients  of 
this  equation  have  like  signs,  and  it  follows  from  Descartes0 
rule  of  signs  that  there  are  no  positive  values  of  j  satis- 
fying this  equation.   Therefore,  the  function  k/e  is  mono- 
tonic  for  j  >  0,  and  a  simple  check  at  j  =  0  and  j  •♦  oo  veri- 
fies that  k/e  monotonically  decreases  to  its  asymptotic  value 

Lim  k,e  =  — — 

It  follows  that  the  critical  value  of  k,e  is  determined  by 
the  mode  corresponding  to  j  ■+  °° .  As  a  result,  the  sign  of 
the  deflection  as  k,e  approaches  the  critical  load  becomes 
indeterminate . 
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The  explanation  of  this  effect  is  that  as  the  number 
of  buckling  waves  increases,  the  bending  stresses  in  the 
faces  will  dominate  the  membrane  stresses.   However,  terms 
due  to  bending  effects  were  neglected  in  the  derivation  of 
equations  (3.44c)  and  (5.1).   Finite  roots  should  result  if 
these  terms  are  retained. 

As  a  result  of  the  preceding  discussion,  it  is  neces- 
sary to  make  a  qualifying  remark  about  the  deflection  curve 
in  Figure  6  for  re   =  0.5,  1.0,  and  2.0.   According  to  this 
discussion,  the  portions  of  the  curves  near  the  buckling 
load  should  be  indeterminate.   The  values  shown  were  actu- 
ally calculated  from  an  infinite  series  truncated  at  the 
term  j  =23. 


CHAPTER  VI 


SUMMARY 


The  equations  governing  the  deflections  of  a  sandwich 
panel  subject  to  a  general  system  of  loads  and  a  thermal 
gradient  were  developed  and  solved  for  the  special  case  of 
uniaxial  loading.   Application  of  a  Levy- type  solution  en- 
abled the  problem  to  be  solved  without  making  simplifying 
assumptions  about  the  boundary  conditions.   The  solution  for 
the  vertical  deflection  was  evaluated  on  a  digital  computer. 

It  was  seen  that  when  bending  stresses  in  the  faces 
are  neglected,  the  response  of  the  panel  is  determined  by 
three  nondimensional  parameters.   The  buckling  loads  were 
seen  to  increase  as  the  aspect  ratio  increased  and  to  de- 
crease as  the  parameter  re  increased.   Reversal  of  deflection 
generally  occurred.   It  was  noted  that  buckling  is  charac- 
terized by  a  large  number  of  half-waves  for  panels  with 
properties  and  dimensions  commonly  used  in  practice.   An 
unusual  effect  of  an  indeterminate  buckling  mode  was  noted 
and  accounted  for  by  the  fact  that  the  bending  stresses 
should  be  included  for  high  buckling  modes. 
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APPENDIX 

The  dimensions  of  the  sandwich  panels  used  by  Ebcioglu 

24 
et  al.    and  the  eight  points  on  the  faces  at  which  the  tem- 
perature was  measured  are  shown  in  Figure  10. 


0.5" 


Figure  10.   Dimensions  and  data  points 
of  the  sandwich  panel 
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The  temperature  gradients  that  were  measured  are 
given  in  Table  4. 

TABLF  4 
EXPERIMENTAL  TEMPERATURE  GRADIENT 


Data 
Point 

Temperature 

Gradient 

(°F) 

Number 

Trial  1 

Trial 

2 

Trial 

3 

Trial  4 

1 

12 

12 

16 

18 

2 

16 

19 

21 

33 

3 

14 

19 

19 

30 

4 

35 

37 

39 

58 

5 

53 

56 

60 

97 

6 

53 

59 

63 

97 

7 

48 

51 

58 

80 

8 

26 

23 

21 

30 
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